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PREFACIO A LA EDICIÓN EN LENGUA ESPAÑOLA 


Con motivo de la decisión de publicar esta obra en lengua 
española desearía dirigirme al lector con algunas palabras. 

Hace ya cerca de veinte años que apareció este libro, En aquel 
entonces la programación lineal se encontraba en el proceso 
de su desarrollo, estábamos apasionados con las cuestiones de su 
teoría, aplicación y divulgación En este pequeño libro traté 
de explicar los conceptos básicos de la asignatura y la esencia 
de algunos de sus métodos de cómputo. 

El libro se agotó rápidamente, continuaron las investigaciones 
sobre el tema; han'aparecido muchas obras en las cuales, bajo 
mi opinión, se explica mejor parte del material. A pesar de eso, 
„hace relativamente poco tiempo, este libro fue reeditado en Japón, 
en la RDA, hi traducido a las lenguas de otros países. 
Ahora la editorial “MIR” está preparando en lengua española 
la edición de la serie completa de libros “Lecciones de divulgación 
de matemáticas” en la que éste salió. La traducción de este libro 
se ha hecho integramente por dicha edición. Observaremos que 
desde aquel tiempo las computadoras se han perfeccionado 
muchísimo у рог eso las referencias а los ejemplos de rapidez 
de resolución de algunos problemas, а} como las dimensiones 
de éstos, deben ser interpretados por el lector como factores que 
caracterizaron el período inicial de la aplicación de la programación 
lineal. Actualmente las dimensiones de los problemas que se tratan 
son mucho mayores y las velocidades de su resolución mucho más 
altas. Las explicaciones que se dan sobre los conceptos y definiciones 
fundamentales de la programación lineal y las bases de los métodos 
de cálculo por el criterio de costo y de tiempo, pueden hoy 
en día abrir al lector el camino a esta interesante asignatura 
de la matemática moderna, a la programación lineal. 

Estaré satisfecho de que el lector encuentre útil para sí este 
pequeño libro y al mismo tiempo le expreso mi simpatía. 


5 de mayo de 1976 A. S. Bársov 
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PREFACIO 


En este libro "se examinan cuestiones de la teoria y de los 
métodos de resolución de algunos problemas de programación 
lineal. El está destinado a un amplio círculo de personas ocupadas 
en el empleo de métodos matemáticos en la organización 
y la planificación de la industria. 

Se estudian los fundamentos de la programación lineal. 
AL hacerlo se presentan tan sólo los datos y demostraciones que 
son necesarios para una exposición elemental de los métodos 
de programación lineal. 

El trabajo se realizó a base de las conferencias dadas por 
el autor en el año 1957 para personas que se dedicaban 
а la resolución de problemas de programación lineal en las 
máquinas computadoras electrónicas. 

L. A. Lyusternik, miembro correspondiente de la А.С. de la 
U. R. S. S, en el año 1959 revisó con atención el material de las 
conferencias, dio una serie de valiosos consejos y contribuyó 
a la edición del presente trabajo. й 

El autor agradece а los'profesores А. A. Liapunov у N. 5. Kra- 
silnikov por su ayuda en la solución de las dificultades que 
surgieron en el proceso de la preparación de este libro. 

El autor está particularmente agradecido al redactor V. D. Ro- 
senko por su minuciosa labor que sirvió considerablemente para 
mejorar el libro. 


А. $. Bársov 
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INTRODUCCIÓN 


La tarea del desarrollo ulterior de las fuerzas productivas, 
de la mejora de la planificación de la industria socialista y el 
aumento de la efectividad económica de las inversiones básicas 
en nuestro país adquiere de año en año cada vez mayor 
importanci 

La diversidad de posibles soluciones técnicas y caminos 
de desarrollo en la industria actual, las interrelaciones entre las 
diferentes ramas de la economía nacional y otros problemas 
económicos hacen que las tareas planteadas anteriormente sean 
excepcionalmente difíciles. 

Para solucionar estos problemas, los métodos matemáticos 
y en particular el método de programación lineal, así como los 
medios técnicos modernos, las máquinas computadoras electrónicas 
pueden prestar una ayuda sustanci 

La teoría de la programación lineal, que su 
décadas, actualmente ha obtenido una amplia utilización práctica 
particularmente en el terreno de la organización y planificación de 
la industria. 

Los primeros trabajos en este sentido fueron los de L. У. Kan- 
torovich, miembro correspondiente de la A.C. de la U.R.S.S.*. 
En esos trabajos se expusieron métodos matemáticos para resolver 
problemas tales como el del aumento de la efectividad del 
transporte, el cálculo de los regímenes óptimos de producción, 
la distribución racional de los materiales industriales, etc. 

Posteriormente fueron creados métodos generales de рго- 
gramación lineal tales como, por ejemplo, el simplex, el combina- 
torio y otros métodos que se emplean eficazmente en la solución 
de diversos problemas para determinar el óptimo en la planifi- 
cación. Dantzing, Charnes y una serie de científicos soviéticos 
y extranjeros se ocuparon de la elaboración de estos métodos. 

La programación lineal abarca métodos de solución de problemas 
de óptimo en los que hay muchas ,variables relacionadas entre 
sí y subordinadas a unas determinadas condiciones de restricción. 
El planteamiento de los problemas de programación lincal se puede 
formular de la siguiente manera: 

Tenemos cierta magnitud (por ejemplo el costo, el tiempo) 


* L, V. Kantorovich hoy es membro efectivo de la 
A. C. de la U.R.S.S. y laureado con el Premio Nobel. 
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que es función lineal de una seric de variables. А su vez, las 
variables tienen que satisíater a las restricciones expresadas en 
forma de un sistema de desigualdades o igualdades lineales. 

Hay que buscar aquellos valores no negativos de las variables 
con los que la magnitud que sea su función adquiera el valor 
menor (el mayor). 

En calidad de ejemplo veamos el problema de transporte. 
Este problema se puede formular de la forma siguiente: 

De т puntos de partida dados, en саба uno de los cuales 
hay a, unidades de cargamento, se deben llevar b, unidades 
de la carga indicada a cada uno de los n puntos de llegada 


(im 1, 2, ‚т; J=1,2 в). 


Se exige planificar el transporte de tal modo que los gastos 
sean mínimos. Designaremos con xy la cantidad de carga trans- 
portada del punto de partida Е al punto de llegada j. En tal 
caso el contenido matemático del problema se reduce a encontrar 
los valores по negativos de xy que satisfagan las ecuaciones 


у= bp 


Ga 


y con los cuales el costo general del transporte 


C= È E, cu 
sea el menor. Aquí су es el costo de transporte de una unidad 
de carga desde el punto de partida i hasta el punto de llegada j. 

En un serie de casos prácticos importantes, el problema se 
plantéa así: se necesita planificar el transporte de cargamentos 
desde m puntos de partida dados hasta n puntos de llegada 
de tal forma que todas las operaciones de transporte se realicen 
en un plazo mínimo. 

Otro ejemplo del empleo de la programación lineal puede ser 
el problema siguiente: 

En muchas fábricas, la producción de diferentes artículos se 
realiza en líneas automáticas. En estos casos pueden surgir diferentes 
problemas referentes а la organización más racional de la producción. 

Supongamos, por ejemplo, que un taller dispone de m máquinas 
para fabricar л artículos distintos. Cada máquina i (i = 1, 2, .... т) 
se caracteriza por et tiempo posible de trabajo mensual 6, 


п 


la norma де tiempo (у para elaborar una unidad del articulo j 
(=1, 2, ~, п) y por el costo су de fabricación de una unidad 
del artículo j en la máquina i. Si al taller зе le ha dado la tarea 
de sacar en el mes próximo una determinada cantidad а, de cada 
uno de los diferentes artículos, surge el problema de organizar 
el trabajo en tal forma que esta tarea se cumpla con los gastos 
mínimos. Si designamos con xy la cantidad de artículos j que 
se fabrican en la máquina i, el problema planteado se reducirá 
a una distribución del trabajo entre las máquinas que satisfaga 
las condiciones 


Aura А tx < bı 


y que lleve el valor total de los costos 


al mínimo posible. 

En los problemas de programación lineal las condiciones que 
se aplican a la zona de valores admisibles de las variables se 
determinan por un sistema de desigualdades o igualdades lineales, 
Al mismo tiempo, la función cuyo valor mínimo (máximo) se busca, 
es también una función lineal de las mismas variables. Esto lo 
subraya la propia denominación de programación lineal. 

Para la determinación de la solución óptima, los métodos 
de programación lineal exigen el examen de varias soluciones. 
Al analizar problemas prácticos, por ejemplo, el problema del 
empleo racional de la maquinaria o de las empresas con determinadas 
condiciones restringentes, el paso de una solución a otra corresponde 
а un estudio sucesivo de diferentes programas de producción. 
De aquí viene el nombre de programación lineal. 

El problema de la programación lineal ез un problema 
de búsqueda de un punto de cierta región en el que la función 
adquiere el valor máximo (mínimo). Por eso surge la pregunta 
natural: ¿Por qué no es suficiente en este caso emplear los 
conocidos métodos clásicos de resolución de los problemas de 
cálculo de valores extremos de las funciones, por ejemplo, el 
método de Lagrange? 

La causa reside en que los métodos clásicos exigen la existencia 
de derivadas parciales de la función en el punto en el que se 
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alcanza el extremo mientras que la función lineal llega a su valor 
extremo en el límite de la región, donde las derivadas parciales 
no existen. 

Eso fue lo que sirvió de motivo para la creación de nuevos 
métodos de resolución de problemas de extremo entre los que se 
encuentra la programación lineal 

La práctica de resolución de problemas de programación lineal 
muestra que cuando es grande cl número de variables, para 
resolver tales problemas es necesario emplear computadoras 
electrónicas. La máquina resuelve en unos dos o tres minutos 
problemas en los que el hombre gastaría hasta una semana, 
Al ser muy grande el múmero de las variables, estos problemas 
pueden ser resueltos solamente por medio de computadoras 
electrónicas. 

Puede servir de ejemplo el cálculo del plan óptimo de transporte 
de arena de construcción a las obras de edificación de la ciudad 
de Moscú. En este problema había 10 puntos de partida y 230 puntos 
de llegada. El plan óptimo de transporte calculado en la computadora 
“Strelá”, en diez días del mes de junio de 1958 dio un 11% 
de economia. 

A continuación se examinan los fundamentos matemáticos y los 
métodos de resolución de algunos problemas de programación 
lineal y, en particular, de los problemas de transporte. 
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CAPÍTULO 1 


ALGUNOS CONCEPTOS Y DEFINICIONES 
DEL ALGEBRA LINEAL 


En este capítulo se exponen los conceptos básicos y las 
definiciones del álgebra lineal de espacio m-dimensional que son 
necesarios para la resolución de problemas de programación lineal. 


$ 1. CONCEPTO DE ESPACIO M-DIMENSIONAL 


Cualquier trio ordenado (a,, аз, ау) de múmeros reales puede 
ser interpretado geométricamente сото un punto en el espacio, 
En concordancia con esta representación geométrica se ha tomado 
en la matemática la definición siguiente: el espacio tridimensional 
es el conjunto de todos los posibles trios ordenados (а, аз, аз) 
de números reales. También se dice que el sistema de números 
(а,, аз, ау) determina el punto М en un espacio tridimensional 
con las coordenadas ay, ау, a о el vector Р con los componentes 
а, аз а; en ese mismo espacio. Estos tres números reales son 
insuficientes para representar ciertos objetos, procesos o estados. 
Por ejemplo, para determinar la posición de un cuerpo sólido 
en el espacio son necesarias seis coordenadas. En el caso de que 
en una zona económica se produzcan determinados artículos 
industriales y agrícolas, como vagones de ferrocarril, automóviles, 
trigo, leche, cerillas y otros, entonces para la característica de esa 
producción industrial y agricola de la гопа se necesitará una 
sucesión ordenada de números reales. Por ejemplo, la tabla 1 
indica que la zona Ne 2 produce anualmente az, toneladas de carbón, 
a toneladas de mineral de hierro, az toneladas de acero, ..„ аз, 
toneladas de trigo. 

Similarmente se definirán, también por medio de sucesiones 
ordenadas de números, por ejemplo: la cantidad de combustible 
de aviación de diferentes clases que se emplea en tal pais, 


* Por analogía el conjunto de todos los posibles 
números reales (а) es un espacio unidimensional cuyo modelo geométrico 
puede ser la recta; el conjunto de todos los posibles pares de números 
reales (ал, аз) es un espacio bidimensional cuyo modelo geométrico puede 
ser la superficie plana. 
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la cantidad de artículos de diferentes clases que se encuentran 
en tal almacén, ete. 


опон % | w | ва 


Estos ejemplos muestran la utilidad del estudio de todos Јов 
posibles conjuntos de sucesiones ordenadas de m números reales 
рага el caso de que т sea cualquier número entero positivo. 
Como se sabe, un sistema ordenado de m números reales 
(а, аз, га Gp ... Gn) se Пата vector m-dimensional. Los números а, 
i = 1,2, .., т, se denominan componentes del vector Р (а, , az, ~., an). 

Los vectores Р (ај, 02, . aa) y Q (bi, Da, ~ bn) se consideran 
iguales si y sólo si coinciden sus componentes que se encuentran 
en los mismos lugares, o sea si a, = Б, para todos los i = 1, 2,.„.т. 

En el caso de que nos interesara conocer la productividad 
total de las diversas clases de producción de dos zonas económicas, 
evidentemente ésta podría calcularse sumando las productividades 
respectivas de estas zonas.* Si la productividad de la zona № | 
en todas las clases de productos se determina por medio del vector 
Р, (а, ала, -~ Gin) y la de la zona № 2, por el vector Р; (а, 022, 
чә Gam), entonces la productividad total de estas dos zonas estará 
caracterizada por el vector 


Оба + аи, аз + а, о Фуа + аза). 


Si la productividad de una zona se determina рог el vector 
P = Р(а,, аз, .. Gn), al incrementarse la productividad de cada 
uno de los productos en k veces, el aumento de la productividad 
de la zona se podrá expresar con el vector Q, el que será 


Q = 9 (Ка), Каз, m, Кал). 


+ Comprenderemos en este caso por productividad 
la cantidad de producto obtenido en un determinado lapso. 
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Las definiciones introducidas son sencillamente una generalización 
de las operaciones ya conocidas con los vectores en un espacio 
tridimensional. Al aplicar el concepto del espacio tridimensional 
а la sucesión de números reales т, obtendremos la importante 
definición siguiente: el conjunto de todos los posibles m-dimensionales 
vectores Р (а, ау, .., am) Con componentes reales se llama espacio 
m-dimensional y se designa Ре". Según la definición, sumar dos 
vectores m-dimensionales consiste en obtener un tercer vector 
R=P+Q con sus componentes iguales a las sumas de los 
respectivos componentes de los vectores sumados. Multiplicar el 
vector P por el número k quiere decir multiplicar cada componente 
por ese número. 

Se dice que el vector Р (а, az, ~., а,) es proporcional al vector 
Оф, ba, m bp) si existe un número k tal que b; = Ка, 
bz = Каз, .., 6, = Ка„ En este caso, P = КО. 

El concepto de la combinación lineal de los vectores es la 
generalización de la noción de proporcionalidad de los vectores. 

El vector P se llama combinación lineal de los vectores Ру, 
Ру, 2, Р,, si existen algunos números reales М, la, .., l, tales 
que Р = ПР, + Р +... + l,P,. En este caso el iésimo componente 
del vector P (i = 1, 2, ... т) será igual a la suma de los productos 
de la multiplicación de cada ¡-ésimo componente de los vectores 
Р, рог su lı, la, .., l, respectivo. 
de vectores Py, Ру, ..., Р, 1, P, se llama linealmente 
dependiente si por lo menos uno de ellos es una combinación 
lineal de los demás vectores. 

Esta definición es equivalente a otra: 

un sistema de vectores se llama linealmente dependiente si 
existen tales números reales kı, kj, ., ky, de los que por 
lo menos uno no sea cero, que satisfagan la igualdad 


КР, + ЮР; +... + КР, =0. 
En el caso contrario, el sistema de vectores se llama linealmente 
independiente. 

Si el vector Py es una combinación lineal de los vectores Pi, 
Ра, ма Ра, se dice que Po se expresa en forma lineal por medio 
del sistema de vectores {Ру}, en el que j = 1, 2, ..‚ п. Está claro 
que si el vector se expresa en forma lineal рог médio de cierto 
subsistema del sistema dado, él se expresará también en forma 
lineal por medio del sistema para lo que es suficiente tomar los 
demás coeficientes del sistema iguales а cero. 
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Generalizando esta terminología, se dice que un sistema 
de vectores Qi, Qi, -~ Q, se expresa en forma lineal por medio 
de un sistema de vectores Ру, Py, ~, P, si cualquier vector Q,, 
i=l, 2, m 5 es una combinación lineal de los vectores del 
sistema {Р} j = 1, 2, .., п. 

Examinemos en el espacio РУ los vectores: 


i (1, O, O, „n 0), 
ia (0, 1,0, ... 0), (0 


1.00, 0, O, ~, 1). 
Estos vectores se llaman vectores — unidad (o vectores — unitarios). 


El sistema de vectores (1) es linealmente independiente, puesto que 
kih + kalz +. + Kaly =0 solamente si k, =0 para todos los 


Po 


0 


I=1,2, m, m. Cualquier vector Р(аз, аз, .., a) del espacio Ре? 
se expresa en forma lineal por medio de los vectores del sistema (1), 
о sea, 

P = ah + asia +. + anin: 


Se puede demostrar que cualquier sistema de vectores de un espacio 
Ре? que esté compuesto por más de т vectores es linealmente 
dependiente. 

De este modo tendremos que і en una superficie salen del 
origen de las coordenadas dos vectores Р, у Р; que по se 
encuentran en una misma dirección, o sea que son linealmente 
independientes, cualquier tercer vector se podrá representar como 
una combinación lineal de estos vectores. 

Por analogía, si en un espacio -tridimensional se dan tres 
vectores qué no se encuentran en шпа misma superficie y que 
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salen del origen de las coordenadas, cualquier vector de este 
espacio se expresa como una combinación lineal de estos vectores. 
En la fig. 1 se ilustra el caso en el que el vector Ро se representa 
por medio de una combinación de los vectores linealmente 
independientes P,, Ра, Py, que tiene la notación siguiente 


1 2 
Po = P, + Pa + Ps 
En un espacio bidimensional, a dos vectores linealmente 


indepándiantes Pi (ага, DIY Palais. exi) les: comréspondo la 
determinante 


ап. Mn | diferente de cero. 
ал an 

El valor absoluto de esta determinante es igual a la superficie del 
paralelógramo construido con los vectores Р, y Р, (fi 
En un espacio tridimensional, tres vectores linealmente indepen 
Р, (аза, аза, аз), P3 (a12, аза, азз) у Рз (ауз, 41, ауз) forman 
un paralelepípedo (fig. 3). En este caso, el valor absoluto de 
la determinante será 


аһ аз аз 
аһ da аз 
ау dy ауу 
diferente de cero e igual al, volumen del paralelepipedo. 


Análogamente, si en un espacio m-dimensional se dan т 
vectores linealmente independientes 


Py (Gr), аз, о, ау, ај), 


entonces, como se demuestra en los cursos de álgebra superior, 
la determinante 


аһ а: im 


an an аы 


айа бв)... Amn 
es diferente de cero. 


2-0307 
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Supongamos que en un espacio m-dimensional están dados n 
vectores cualesquiera 


CA ар ан) 


Formemos соп los componentes de estos vectores la matriz” 
de un orden-de (т х п) 


„з тез 
P, Р; Р, Р, 
аа ап ау а, 
Ia а - ву в 

ч @ 
РИН а 
ба аа. а 


Hablando en general, las columnas de esta matriz vistas como 
vectores m-dimensionales pueden ser linealmente dependientes. 
El máximo número de columnas linealmente independientes de 
la matriz (2) se llama rango de esta matriz*". Dicho de otra 
forma, el rango de la matriz (2) es igual al máximo número 
de vectores linealmente independientes Ру, cuyos componentes 
forman sus columnas. Cualquier sistema linealmente independiente 
máximo de vectores del espacio P™ se Ната base de ese espacio 
vectorial. 


+7 Aqui y а continuación colocaremos los componentes 
de los vectores Ру en forma de columnas de las matrices. 
El número máximo de filas linealmente indepen- 
нез de la matriz siempre es igual al número máximo de columnas 
linealmente independientes. 
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Supongamos que los vectores Py, Pz, .., Ру-1, Ру Ру, => Ра 
forman la base del espacio РЕЧ la que llamaremos A-base. 
Entonces cualquier vector P de ese espacio puede ser representado 
en la base dada en forma de una combinación lineal de vectores 
Py j= I, 2, ... т. Además, esta forma de representación tendrá 
una sola aceptación. Todas las bases de un espacio vectorial se 
componen del mismo número de vectores. 

“Tomemos un vector cualquiera Q que no pertenezca a la base 
Py, Ру, ~, Pa. Entonces, si el vector Q no es igual a cero, 
en la combinación lineal 
Ф = о.Р, + аР +... + Y Pj-1 + аР + ау Ра +. GP 
habrá рог lo menos un coeficiente diferente de cero. 

Supongamos, por ejemplo, que ај 0, entonces el vector Р, 
de la A-base se puede representar de la siguiente manera 


Excluimos de la A-base el vector Ру y añadimos a los vectores 
que quedan el vector Q. El sistema de vectores Ру, Pa, .., Py 45 
О, Ру», ~» Ра de nuevo es base. Realmente, teniendo en cuenta 
que la descomposición de cualquier vector entre los vectores 
de la base se efectúa siempre de una sola forma, después de la 
exclusión del vector Ру de la A-base, el vector Q ya no podrá 
ser representado como una combinación lineal de los demás 
vectores. Eso quiere decir que el sistema de vectores Pi, Pa, .., 
Ру, О, Рац, => Ра 65 un sistema lincalmente independiente, 
о sea, que forma una Базе. 

La B-base, que se ha obtenido de la А-ђазе por medio 
de la sustitución de cierto vector por un vector no perteneciente 
a la А-Базе, la llamaremos base соп una sustitución con respecto 
a la A-base. 

Supongamos que en un espacio m-dimensional están dados 
л+ 1 vectores Po, Ра, Pz, -~ Ру, - Ра. los que en la base 
Qi, Qz, .., Qu tienen la descomposición 


Po =b,Q, + 5,0; ++ ВА Фа 
Р, =а110, + 21,0; +...+ ami да 
Ра = а10, + 2,10; +... ад 


B 


41,0, + а;/0; +...+ а, 


Р, = 01,0, + 1,0 а 


Р, 
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Formemos una matriz con los componentes de la descomposición (3) 


A Б 
р. зид ¿do б эз ы 
A A 
E (4) 
ы а а а, 


CARA A ж 


Admitamos ahora que los vectores qı, Фу, ~» Фа forman una 
nueva base P™, Supongamos que estos vectores tienen la siguiente 
descomposición en la base Qi, Qa, ... Qu: 

qı =411Q1 + 4110; + ~ + да Qu, 

9 = 420. +4220: + „+ Gm: Om, 


9. = 4.01 + ФО: + ~ + Ion а 


Сото qı, 92, = Фа forman una base, la determinante 


Qu Qz o Фе 
Gar 92 o. Фа 
a=l: : 6 


Ami na зс Amm 


será diferente de cero. 

Pasemos de la base Qi, Qz, m Qu а la base qı, 92) ~ Че 
Entonces, los vectores Po, Ра, Ра, .., Р, tendrán сп la base 
di> 4, — Чл Unos coeficientes de descomposición que, en general, 
serán diferentes de los coeficientes de descomposición en la base 
Qis ©, ~ Qu: 


Po = biqi + Баба + ~. + Бф, 
Ру = ауф + 43/92 +... + ара Û = 1, 2, „~ 1). 


En este caso, la matriz de la descomposición ‘tiene la siguiente 


21 


notación 
Po Р, P, 
ы ац з din 
ГА бб a 
: : © 
K а % а 


6 чене ФЬ, 


Сото se sabe, los elementos de la matriz (6) se determinan por 
medio de los elementos de la matriz (4) y de la determinante (5) 
con las fórmulas siguientes: 


Я ‚ Aj 
لك‎ a, [0] 
en las que 

Qu ят 
Аў = | 91 Фа 

Ami Amm 

Qu im 
Af = | 41 qam 

Gm ds 

(=1, 2, .. 


A continuación, al mismo tiempo de denominar estas determinantes 
con Д, Af, Af, emplearemos a veces las designaciones siguientes: 
А = (19 — qm); АР = (а ~ 
А = (ф ~ Puss - 4). 
Examinemos un ejemplo de determinación de la descomposición 
de un sistema dado al pasar a una base nueva. 


Supongamos que en la base P,, Ра, Pa, los vectores Pò; Ps, 
Pa, Pa, Pa, Ps, Ре tienen unas descomposiciones que se determinan 


22 


рог la matriz 


5 
3010 –-:1 
2 


-2 
KEFTE 4 
Puesto que la determinante compuesta de los componentes de los 
vectores Pa, Ps, Ре 
-3 4 -1 
А = (Ра, Ps, P) =| -1 1 -2|=-18 
2d 4 


es diferente de cero, los vectores Pa, Ра. Pg formarán una base. 
Encontremos la descomposición de los vectores Po, P,, Pa, Py, 


Pa, Ps, Ps еп la base Pa, Ps, Pg. Empleando las fórmulas (7), 
obtendremos 


Po Р, Pa Ру 
(PoPsPeh (Р,Р;Ре), (P2PsPo) (PsPsPo), 
(P.PoPo), (PPPs) (РаР2Р5), (PaP3Po), 
(PaPsPo), (РРР), (PaPsP2), (PaPsP3), 


Pa Ps Pe 
(PgPsPe), (PsPsPo) (PsPsPe) 
(РАРАРе), (PaPsPo), (Р„РЬР„) 
(Р.Р,Р,), (PaPsPs) (PaPsPo) 
Realizando las transformaciones correspondientes, obtendremos 


la matriz que representa la descomposición de estos vectores 
en la base Pa, Ps, Ре: 


Р, P, Ра Ру P, 


= 1 0 0 
18 18 

| EE ê i 
18 18 

-i L 0 e 
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Como es sabido, se llama producto escalar de dos vectore: 
P y Q ha expresión definida por la fórmula Я 


(Р.О) = a,b, + а,Ь, +... + abi +... + а,Ь, 


donde a, b son los componentes de los vectores Р у Q 
Los vectores P y Q se llaman ortogonales si su producto escala 
es igual a cero. Ya que зе puede considerar a una matriz de ш 
orden de mx п como un vector (т х n)-dimensional, la Шата 
remos producto escalar de las matrices 


aa ж =й 
ан аз 
A= $ 
C و ق‎ аы 
у л 
bu ba 
Ба ба 
B =| : 


bm ba 


a la suma algebraica de todos los elementos de la matriz 
аһ bu an ba а. bie 
da ба an ba +з аһ Фа 


аа. De гаа Va we Дш б 


y la designaremos (А. В). 
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$ 2. HIPERPLANO Y SEMIESPACIO 


Por la Geometría Analítica sabemos que a la ecuación lineal 
Ајх + 4х: + Духу = С (8) 
en un espacio tridimensional le corresponde el plano normal 


al vector А(Аз, Az, A 
Reduzcamos la ecuación (8) a la forma 


EE + 22. 1. (9) 
а aa 
El plano que corresponde a esta ecuación, corta en los ejes 
de las coordenadas los segmentos ау, аз, ay (fig. 4). 


O-origen de las coordenados 


а mes 


La ecuación de una superficie plana en un espacio tridimen- 
sional puede ser, además, representada en forma vectorial 
(А.Х) = h, donde А“ es un vector de longitud igual a uno 
(unitario), normal al plano; X, el vector corriente que une el 
origen de las coordenadas con cualquier punto perteneciente 
al plano. Aquí (А" - X) es el producto escalar de la multiplicación 
de los vectores A” y X que tiene un valor igual a la proyección 
h'del vector X en la dirección determinada por el vector A”. 
El valor de h es igual a la distancia entre el origen de las 
coordenadas y el plano. El plano pasará por el origen de las 
coordenadas cuando л = 0. 

Cualquier vector Х que una el origen de las coordenadas соп 
un punto del plano tendrá la misma proyección h en la dirección 
A” (fig. 3. 
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Análogamente llamaremos hiperplano, o sencillamente plano 
en el espacio m-dimensional, al conjunto de todos los puntos 
ба» ха, => Xm) que satisfagan la ecuación — * 

Ауху + Ауху фи. + А„х„ = С. (8) 
Convendremos en decir que este hiperplano es normal al ‘vector 
A(Ay; Az, ~, An). Diremos también que la ecuación 


A ПИР س‎ (у 
а а Гад 


corresponde al hiperplano que corta, en los ejes de coordenadas, 
segmentos de longitud igual а а, аз, = ал. 


2 


Recta Р 


Ecuación de la recto Хь " 
humos 


мо ASA 


Ecuación del plano 
-51t r= 31, =5 


6 sig 7 


Por último, la ecuación vectorial (A*-X)=h define en el 
espacio Р"? un hiperplano normal al vector unitario A” que se 
encuentra a una distancia h del origen de las coordenadas. 

Una recta en un plano divide a éste en dos partes, cada una 
de las cuales se llama semiplano. En la fig. 6 se ve que una recta 
divide un plano en dos semiplanos. Añadamos que la proyección 
Осу de los vectores X,, cuyos extremos se encuentran en uno 
de los semiplanos, son menores que A = Ос, y la de los vectores 
X,, cuyos extremos se encuentran en el otro semiplano, son 
mayores que h. 

Un plano en un espacio tridimensional también divide todo 
el espacio en dos partes, cada una de las cuales se denomina 
semiespacio. Й 

Por analogía, diremos que un hiperplano en un espacio 
m-dimensional divide este espacio en dos partes cada шпа de las 
cuales se denomina semiespacio. 
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Supongamos que un hiperplano en с! espacio P™ se expresa 
por medio de la ecuación (A*-X) =h. Entonces, para los puntos 
М de uno de los semiespacios, las proyecciones que representan 
sus vectores X en la dirección А° serán menores que h, y para 
los puntos del otro semiespacio, mayores que A. Así diremos que 
uno de los semiespacios es el conjunto de vectores X para los 
que se cumple la desigualdad (A*-X) < h, y para los vectores del 
otro semiespacio, la (A*-X) > А. El propio hiperplano (A°- X) = h 
puede ser agregado a uno de los semiespacios. Entonces to 
el conjunto de puntos del espacio m-dimensional se divid 
en dos partes: la de los puntos para los que (A*:X)<h 
y la de los puntos para los que (А.Х) > h (о bien, en otro caso 
(A°-X) < h y (A°-X) > h). 

Ejemplo 1. Dado el semiespacio — Sx, + 4x — 3x, < 5 deter- 
minar si le pertenece el punto (0, 0, 0). 

Para obtener la respuesta es suficiente colocar los valores 
х =0, ха =0, x3=0 en la desigualdad. Tendremos —5.0+ 
+4-0—3-0<5, de lo que se deduce que el punto (0, 0, 0) 
pertenece realmente al semiespacio —Sx, + 4x3 — 3x3 < 5. 

El plano —5ху + 4x — 3x3 = 5 es perpendicular al vector 
А(—5, 4, —3) (Mg. 7). 

Ejemplo 2. Determinar si el punto de un espacio 9-dimensional 
xı =0, x =4, x3=3 х= —7, х = 0, xe =0, ху =9, ху = 1, 
ху = 0 pertenece al semiespacio 4х, + Sx — 7ху + xs — 2х+ 
+ 12x3 — 3ху < 11. Colocando las coordenadas de este punto en la 
desigualdad, obtendremos 4-0 + 4.5 – 3.7 = 7.0+0:1-2.0+ 
+9۰12 +0۰1 — 3-0 = 107 > 11. En consecuencia este punto по 
pertenece al semiespacio dado. 


$ 3. POLIEDROS CONVEXOS 


Se llama cuerpo convexo a un cuerpo que contenga junto con 
dos cualesquiera de sus puntos, todo el segmento rectilineo. que los 
une. 

Ejemplos de cuerpos convexos son un circulo, una esfera, 
un cubo, un ángulo formado por dos líneas que salen de un 
punto (fig. 8), ete: 

Supongamos que los puntos x e у son comunes para los 
cuerpos convexos А y В (fig. 9). Entonces х e y pertenecen al 
cuerpo А y por eso el segmento que une los puntos x e y 
también pertenece a A. Por analogía, este mismo segmento 
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pertenece también al cuerpo B. Por consiguiente él también 
pertenece a la parte común de los cuerpos A y В. Esto significa 
que la parte comin о la intersección de los cuerpos сопрехоз es 
un cuerpo convexo. 

En una superficie plana tomemos un polígono que se encuentre 
siempre a un lado de cada recta que forma este polígono. Este 
polígono, como se УБ en la бр. 10,a, es convexo. Realmente, 
dos puntos cualesquiera x e y de este polígono le pertenecen 
junto con el segmento que les une. Al contrario, el polígono 
representado en la fig. 10, b no se encuentra por entero en un lado 
de cada una de las rectas que forman este polígono. Un polígono 
tal no es convexo. 


ФА 


“Tomemos un punto cualquiera М que no pertenezca а un 
polígono convexo (fig. 11, a). Siempre se puede trazar una recta PQ 
tal que el punto М y el polígono queden a diferentes lados 
de ella. Se puede construir para un polígono convexo un conjunto 
de rectas tales que cada una de ellas tenga por lo menos un punto 
común соп el poligono y que todo el polígono quede а un lado 
de cada una de ellas. Estas rectas se llaman soportes. 

En el jemplo de la fig. 11,b las rectas AD, CB, DB y FN 
son soportes. La recta soporte puede tener una parte común 
con el polígono convexo. Esta puede ser un punto o un segmento. 

En un espacio tridimensional, un cuerpo limitado por planos 
que se encuentra por entero a un lado de cada plano que 
contiene una de sus caras, es convexo y se llama poliedro convexo. 
Ejemplos de tales cuerpos son un diamante poliédrico, un prisma, 
etc. El plano que tiene por lo menos un punto común соп el 
poliedro y le deja a todo él a un lado de si se llama plano 
soporte de un poliedro convexo. El plano soporte puede tener con 
el poliedro una parte común que puede ser un punto (llamado 
vértice del poliedro), un segmento (llamado arista), o un polígono 
(llamado cara). Está claro que por cada vértice o arista del 


тыз 
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poliedro se puede trazar una cantidad infinita de planos soportes, 
mientras que por сада cara puede pasar sólo uno. 

Anteriormente nos convencimos de que la parte común 
de varios cuerpos convexos es un cuerpo convexo, Por eso, 
la parte comin de varios poliedros convexos es también un cuerpo 
convexo. Ya que el plano es un cuerpo convexo, la intersección 
de un poliedro con un plano también es un cuerpo convexo que 


a A 
аы 


и 


puede ser un punto, un segmento о un polígono convexo. 
Las propiedades de los cuerpos convexos de los espacios multi- 
dimensionales se pueden examinar análogamente a las propiedades 
de Поз cuerpos convexos del espacio tridimensional. Algunas 
de estas propiedades se estudiarán en los $3 4 y 5. 


$ 4. SISTEMA DE DESIGUALDADES LINEALES 


Supongamos que en un espacio bidimensional se dan л 
desigualdades de la forma 
аах, + аах, S by @ф=1,2,.,‚л)%® (10) 
Cualquier desi јад de la forma ах, + ах; > bı 
después de multiplicar E E 
forma (10). 


29 


Cada una de estas desigualdades determina uno de los dos 
semiplanos que tienen por línea límite-la recta apx; + ахз = by... 
La recta-límite ах, + a,x = b; es normal al vector A, (ац, 12). 


Ecuación de: recto 
25936 


ne 


Llamaremos solución del sistema dado a cualquier par de 
números (ху, ха) que satisfaga todas las desigualdades del sistema 
(10). En otras palabras, cualquier punto del plano (xi, ха) cuyas 
coordenadas satisfagan el sistema (10), es una solución. 


Veamos algunos ejemplos: 
1. La desigualdad 


+ <1 û 2x +3 <6 


determina un semiplano (fig. 12, a). A esta desigualdad le satisface 
cualquier punto que se encuentre en la parte rayada del plano, 
La recta-límite se expresa con la ecuación 2х, + 3x7 =6 y es 
normal al vector А (2, 3). 

2. Dos desigualdades 


2x +3x 66 y —х +х,<2 
determinan una parte del plano, como se muestra en la fig. 12, b. 
3. A las tres desigualdades 
2x, +3x, < 6, 
-x+ x1<2 
=x – а 53 
les satisface el conjunto de puntos del plano que forman el 


triángulo de las soluciones AFB (fig. 12, с). 
4. A las cuatro desigualdades 


2х, + Зх, <6, () 
=х +, 52, (2) 
= = 3x 63, o 

х, <4 4 


les corresponde el conjunto de puntos del plano que forma 
el polígono de las soluciones ABCD (fig. 12,2). 
5. А las siete desiguladades 


2x1 + 3x <6, (0 
ови 52, (2) 
-=x = 3х: <3, (3) 
2x1 <3, 4 
= <3 (5) 
3x1 + 7x < 21, (6) 


ж-— За <3 о 
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„les corresponde el mismo conjunto de puntos que en el ejemplo 
Las desigualdades (5), (6) y (7) pueden ser excluidas sin alteración 
del conjunto de soluciones. En este caso, las desigualdades (5) y (6) 
definen las rectas-límite que no tienen puntos comunes,con el 
polígono ABCD. La recta (7) tiene un punto común con 
el poligono y es soporte (fig. 12, e). 

6. El sistema de desigualdades 


2x, +3x <6. aj 
x +X <2, о) 
=x = За 43, (6) 
2x, <3, (4 
—3x, = 2× < —12, (5) 


no tiene ni una solución. Desde el punto de vista de la Geometria, 
esto quiere decir que no existe пі un punto cuyas coordenadas 
satisfagan todas las desigualdades (fig. 12, /). 

El estudio de estos ejemplos поз lleva a las siguientes 
conclusiones: 

1. Un sistema de desigualdades con dos variables puede ser 
compatible. Entonces existe por lo menos un punto del plano que 
pertenece a todos los semiplanos determinados por el sistema dado 
de desigualdades. El conjunto de todos los puntos tales puede ser 
un semiplano, un poligono limitado o ilimitado", una recta 
о su segmento y, por бп, un punto. El conjunto de puntos que 
satisfacen el sistema de desigualdades es un cuerpo convexo, 

2. Un sistema de desigualdades puedo ser incompatible. En este 
caso, no existe ni un punto del plano que satisfaga simultáneamente 
todas las desigualdades del sistema. 

En un espacio tridimensional, un sistema de n desigualdades 
puede escribirse, sin limitarle la generalidad, en la forma 


апх + аах + axs <b, (81,2, п}. an 


Como ya sabemos, cada una de las desigualdades (11) define 
un semiespacio con el plano limite 


аах + аах) + dX = Ы. 


% Pueden haber sistemas de desigualdades а los que 
correspondan figuras (poliedros) no limitados por sus aristas y caras en 
alguna de las direcciones. Estos sisiemas de desigualdades también 
responden a esta propiedad. 
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El conjunto de desigualdades (11) puede ser compatible. En este 
caso, existe cierto conjunto de puntos del espacio tridimensional 
que satisfacen el sistema de desigualdades. Este conjunto es un 
conjunto convexo y puede representarse en forma de un semiespacio, 
un poliedro, un plano, un polígono, una recta о, por último, 
un punto. 

En el caso de la compatibilidad, entre las desigualdades del 
sistema puede haber desigualdades “excesivas”, que al ser eliminadas 
по modifican el conjunto de las soluciones del sistema de 
desigualdades dado. Estas desigualdades “excesivas” pueden ser 
de dos clases, Las de la primera clase son aquellas cuyos planos- 
límite no tienen intersección con el conjunto de todas las soluciones 
del sistema. Las de la segunda clase son las desigualdades cuyos 
planoselimite son soportes рага el conjunto de soluciones. 

Puede resultar que en un espacio de tres dimensiones no haya 
ningún punto que satisfaga simultáneamente todas las desigualdades, 
En este caso, el sistema de desigualdades se llama incompatible. 

Supongamos que еп un espacio m-dimensional'se,da el sistema 
de desigualdades 


аах + арх; Rx Sb (=1, 2 „0 (12) 


Por analogía con el espacio tridimensional diremos que cada 
una de las desigualdades (12) define en el espacio m-dimensional 
un semiespacio соп el plano-límite 


ваха + арха +. + а„х„ = br 


:Si en el espacio т-ателејола] existe por lo menos un punto 
“Мба, Xa, са Хе) que satisface simultáneamente todas las 
desigualdades, el sistema (12) se llama compatible, Al conjunto 
de todos los puntos tales lo llamaremos “el poliedro de las 
soluciones”, 

Supongamos que en un espacio mdimensional se dan dos 
puntos М'(х\, Хр, = XA y M(x, Xi, «a х1). El conjunto 
de puntos M (ха, ха, - Xa) cuyas coordenadas satisfacen las 
condiciones 


х= xj + (x — XD, 
Ху = x3 + (x3 — x), 


ха + (а x) 
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al variar el parámetro г de O a 1, se Пата segmento que une los 
puntos M’ y M”. El “poliedro de las soluciones”, siendo la 
intersección de los semiespacios, es un conjunto convexo. Eso quiere 
decir que junto con los puntos М” у M” también le pertenecen 
todos los puntos del segmento que los une. 

Las desigualdades que pueden ser eliminadas del sistema (12), 
sin modificar el conjunto de sus soluciones, se llaman dependientes 
o “excesivas”. Al eliminar consecutivamente del sistema dado una 
desigualdad de esta clase tras otra, obtendremos un subsistema 
de desigualdades con un conjunto de soluciones que coincidirá 
con el conjunto de soluciones del sistema inicial. 

La eliminación de las desigualdades “excesivas” es un proceso 
muy complicado y trabajoso. Una de las particularidades dé los 
métodos de programación lineal consiste en que para calcular 
el valor mínimo (máximo) de la función lineal en el poliedro 
по es imprescindible eliminar las desigualdades “excesivas”. 
Si al sistema de desigualdades (12) no le satisface ningún punto 
del espacio m-dimensional, este sistema se llama incompatible. 


$ 5. EL VALOR MÍNIMO Y MÁXIMO 
DE LA FORMA LINEAL EN EL POLIEDRO 


Examinemos un sistema compatible de ecuaciones lineales con 
dos variables. Consideremos que ya hemos eliminado todas las 
desigualdades “excesivas” de la 1-а y de la 2-a clase, y asi hemos 
seleccionado el polígono de las soluciones en su “forma pura” 
(fig. 13). Supongamos que, además, se da la función lineal 
de dos variables 


J = суху + саха. 


Encontremos en el conjunto de puntos (xı, ха) del poligono 
de las soluciones tales puntos que lleven» la función lineal 
f = суху + cx, al valor mínimo y máximo. Examinemos el conjunto 
de todos los puntos (xı, ха) del plano en cada uno de los 
cuales la función / = сух, + cx toma un valor fijo (f = f.) 
El conjunto de tales puntos es la recta суху + сах = f ,. Esta recta, 
como se ha señalado en el parágrafo anterior, es normal al vector 
С(с,, сз) que sale del origen de las coordenadas. Tracemos una 
recta F (fig. 13) normal al vector C y desplacémosla paralelamente 
a sí misma en la dirección positiva del vector C. Supongamos que 


1 3.0307 
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en su desplazamiento, la recta Ј se encuentra por primera vez 
con el poligono en el vértice A. En esta posición F' la recta Е 
se hace soporte. Al continuar el desplazamiento en esa misma 
dirección, la recta F pasará por el vértice В y se hará también 
recta-soporte. Puesto que el sentido del vector С(сү, сз) es la 
dirección del máximo incremento de la función lineal f = суху + 
+ саха, entonces entre todos los valores que toma la función 
lineal / en el poligono de las soluciones, esta función tomará 
en la recta-soporte Р' su valor mínimo y en la recta-soporte Е", 
su valor máximo. 


вео ем 


Así que los valores mínimo y máximo de la función lineal 
f = суху + суху en el polígono de las soluciones se alcanzarán 
en los puntos de intersección de este poligono con las rectas- 
soportes normales al vector C(c,, су). La intersección de la recta- 
soporte con el polígono de las soluciones puede ser un solo 
punto” (un vértice del poligono) o un conjunto innumerable 
de puntos (en este caso el conjunto es un lado del polígono). 

En la fig. 14 se presenta el caso en el que la función lineal f 
alcanza el valor mínimo en cada uno de los puntos del segmento 
Ау, Аз, mientras que el valor máximo, en los puntos del polígono 
que están infinitamente alejados. 


# Puede resultar que este punto se encuentre en cl 
infinito, 
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Por analogía, una función lineal de tres variables f = сул, + 
+ саха + суху toma un valor constante en un plano normal 
al vector С(сү, cz, сз). El sentido de С es Іа dirección del máximo 
incremento de la función / (fig. 15). Los valores máximo y minimo 
de esta función en el poliedro de las soluciones también se alcanzan 
en los puntos de intersección de este poliedro con los planos-soportes 
normales al vector С (c1, са, су); aquí la función f alcanza en uno de 
los planos-soportes el valor mínimo y en otro, el máximo, La inter- 
sección de un poliedro con un plano-soporte puede ser un punto 


А 


Pohodro de los 
soluciones 


тын 


(un vértice del poliedro) o un conjunto innumerable de puntos 
(este conjunto será un arista o un lado del poliedro). 

Por ejemplo, en la fig. 16 aparece el caso en el que / alcanza 
su valor mínimo en cada uno de los puntos del lado ABQP 
y su valor máximo, en el punto D. 
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La generalización del concepto de las funciones lineales de dos 
о tres variables es la función f = cx + сх: +. + сух, de n 
variables reales ж, ха, ~ X llamada forma lineal. Aquí са, 
С, m С, SON números reales. Fijando valores de la forma lineal, 
о sea, Јел; /=/; mi f= f, se definen en un espacio 
n-dimensional los hiperplanos 


Гу = сух + суху +. + сух, 
fa = суң + суху +. + сох, 


Sa = аха + саху ++ сох, 


que son normales al vector С(с,, са, ..‚ с,). Designemos con 
Fm el valor mínimo de la forma 


S = соч + аха + + охь 


en el poliedro de las soluciones, y соп fm, el valor máximo. 
Puesto que el vector C determina la dirección del incremento 
máximo de la forma lineal f, entonces cuando las funciones 
f’ < f min Y f’ > f ыз» los hiperplanos correspondientes f’ = сүх + 
+ саха + u. + сух, no tienen intersección con el poliedro de las 
soluciones. Por otra parte, cada hiperplano f” = сух, + саха + 
+... + саху, рага el que fmin < f” < fas, tiene puntos comunes 
соп el poliedro de las soluciones. El conjunto de puntos М (x,, 
Xa, =» Xp) en los que la forma lineal alcanza el valor mínimo 
es la intersección del poliedro de las soluciones con el hiperplano 
soporte сүх + суху +. + их, = Гым, normal al vector С(а, 
са, « со). Análogamente, f alcanza el valor máximo en los puntos 
de intersección del poliedro con el hiperplano-soporte суху + 
+ суху +. + Ca Xp = f qa, que es también- normal al vector С. 

La intersección del poliedro de las soluciones соп un 
hiperplano-soporte será un vértice, una arista o un “lado” del 
poliedro. 

Así que el valor óptimo de la forma lineal en el poliedro 
de las soluciones se alcanza en los puntos entre los cuales siempre 
se encuentra aunque sea uno de los vértices. Por eso para 
calcular la solución óptima es suficiente encontrar el vértice del 
poliedro en el que la forma lineal alcanza el valor mínimo 
(máximo). 
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$ 6. REDUCCIÓN DE DESIGUALDADES A IGUALDADES 
AL RESOLVER PROBLEMAS DE PROGRAMACIÓN LINEAL 


Se da un sistema de m desigualdades lineales con n variables: 


аах + азху + ах, Ы, 
ох + азаха + + ахь ba, 


(2) 


Amy Xy + анха А + амаи Вие 


que determina en un espacio n-dimensional un poliedro de las 
soluciones. Simultáneamente a este sistema de desigualdades veamos 
un sistema de m ecuaciones algebraicas lineales con n + m variables: 


аах фаха + ах р = bi, 
анх + йззХх; +. + Gaga + йу = ba, 


(12) 


Amr Xy + йа + + na + а = Ви» 


Demostremos que a cualquier solución xf, x$, ... x? del sistema 
de desigualdades (12) le corresponde una determinada solución 
Xf, х}, <s XB: Hf, H3» =» Ha del sistema de ecuaciones algebraicas 
lineales (12), teniendo en cuenta que las variables complementarias 
(variables de holgura) satisfacen la condición de que p? > 0, и" > 0, 
~, Из > 0. Realmente, si x, x3, ~, ха es la solución del sistema 
(12), tienen lugar las desigualdades 


AA+ <, 
ах] + ах} + ~ фаха & bz, 


axî + арж ++ аа < bn, 


аах + ах +. + а„ху&Ь„. 
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Designemos 


HÎ = bı — (ах + азах) +. + ах), 
Hî = ba — (аах + 0393 +. + аы), 


(12) 


Entonces, en primer lugar, рӯ > 0, из > 0, .., и» > 0; en segundo 
lugar, como se deduce de (127) el sistema de números xî, x, ~ 
ХИ Hî, Há, - HA ез la solución del sistema de ecuaciones 
algebraicas lineales (12). 

Demostremos, a la inversa, que a cualquier solución xî, xã, 
ха из „ Wa del sistema de ecuaciones (12), que satisfaga 
la condición p? > 0, рі > 0, ..., и > 0, le corresponde una solución 
determinada del sistema de desigualdades (12). En efecto, como 
el sistema de números xî, x 5 Mis р}. =, ир ез la solución 
del sistema (12), podemos escribir: 


Ax + аах} + аи, 


A Xi а: +. + й„х Jl Б, 


De acuerdo a las suposiciones hechas los números pî, Из, 
son no negativos, entonces se cumplen las desigualdades: 


TN 


ах] + ах ++ анху < ба, 
o sea, el sistema de números xî, x, .., x es la solución del 
sistema de desigualdades (12). 

De este modo se ha establecido la existencia de una relación 
“recíproca entre el conjunto de todas las soluciones xf, ха, wa x? 
del sistema (12) y el conjunto de las soluciones xî, х}, mn ~ 
lf, H3, n Hé del sistema (12) en los que se encuentran los mismos 
valores. Observemos que las variables complementarias p$, p3, <= 
из son по negativas. O sea, el problema de la resolución de un 
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sistema de desigualdades lineales (12) se reduce a la resolución 
de un sistema correspondiente de ecuaciones lineales (12). 

En los problemas de programación lineal поз interesarán 
soluciones de sistemas de desigualdades que satisfagan la condición 
xı > Û, ха > 0, ~. х, > 0. Estas soluciones se Патап no negativas, 
Por eso, como la resolución de un sistema de desigualdades se 
reduce a la resolución de un sistema correspondiente de ecuaciones 
lineales algebraicas, una de las tareas más importantes de la 
programación lincal es el cálculo de las soluciones no negativas 
de un sistema de ecuaciones lineales. 

Ejemplo. Calcular la solución no: negativa del sistema de 
desigualdades: 


2х + 3x < 6, 
-x1+x<2 
=x; – Зх < 3. 


Introduciendo las variables complementarias p, > 0, pz > 0, ру > 0, 
obtendremos el sistema de ecuaciones: 


2x1 + 3x2 + и =6 
x + ox + =2 
=x, = а +m=3. 


A cualquier solución no negativa de este sistema de ecuaciones 
le corresponde una determinada solución no negativa del sistema 
inicial de desigualdades y viceversa. Por ejemplo, a la solución 
no negativa х = 1, xı = 1, p= 1, на = 2, py = 7 del sistema 
de ecuaciones le corresponde la solución no negativa x; = 1, x = | 
del sistema de desigualdades inicial. 

Observemos que si el sistema de desigualdades es incompatible 
en la región de las soluciones no negativas, el sistema correspondiente 
de ecuaciones lineales no tiene ni una solución no negativa. 

Es sabido que el valor mínimo (máximo) de la forma lineal, 
en un poliedro determinado por un sistema de desigualdades, 
se alcanza en un determinado vértice del poliedro de las soluciones. 
Se puede demostrar que a cada vértice del poliedro de las soluciones 
le satisface una solución no negativa del sistema correspondiente 
de ecuaciones algebraicas lineales en el que por lo menos una 
variable complementaria es igual a cero. Las variables comple- 
mentarias p; > 0, p3 > 0, ..., ia > 0 que se introducen en el sistema 
de desigualdades se pueden explicar en un sentido geométrico. 

Examinemos el poliedro de las soluciones del sistema de 


desigualdades: 
аха + ах: + ах < by 


беху + 4,32 + ... + ах, S ба 


Supongamos que existe una de las soluciones no negativas xf, 
X$, чч x2 del sistema de desigualdades, A esta solución le corresponde 
el punto M (xî, x$, ~. x?) que pertenece al poliedro de las 
soluciones. 

Supongamos que con h se denomina la distancia entre el punto 
М y el і-ёзіто hiperplano. Los valores de las variables comple- 
mentarias 


ИЙ = bi — (axî + ах} +. + ах), 


= ba (ахї + ах +. + ад 
serán proporcionales a las distancias entre el punto М (xî, x3, ~ 
х) y los hiperplanos-límit 

ах + ах + ах, = 


Amy Xy + а, +. + ах, = ba. 
Realmente, según la definición, la distancia entre el punto М 
y el ésimo hiperplano es igual а 
hm b, ~ (auxi + apx ++ 
| VETA +- + ah 


@=1,2,3, n т), 


de lo que se deduce que: 
мећу +а + ~ + ab, 
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CAPÍTULO П 


RESOLUCIÓN DEL PROBLEMA GENERAL 
DE LA PROGRAMACIÓN LINEAL 


La programación lincal es una disciplina dentro de las 
matemáticas en la que se estudian los métodos de resolución del 
siguiente problema: se da un sistema de m ecuaciones algebraicas 
lineales que contienen п variables: 


аң HG +. "ај + X= у, 
ах + арха + Ф аууху+ + аз„х‚ = Ва, 


анх + азха +. + yx] ++ ах, = bı 


(13) 


+ аах +. + ау А А дахи = Ва 


y la forma lincal f de estas mismas variables 
S = суху + саха + рр + Can (14) 


Hay que encontrar entre las soluciones del sistema (13) una 
solución no negativa en la que la forma lineal f tome su valor 
mínimo”. А esta solución la llamaremos la solución óptima 
de este problema. 

En сме capítulo se estudia un método de cálculo de una 
de las soluciones по negativas del sistema (13) y uno de los 
algorítmos de resolución del problema general de la programación 
lineal. 


$ 7. TRANSFORMACIONES IDÉNTICAS DE 
UN SISTEMA DE ECUACIONES ALGEBRAICAS LINEALES 


Se da un sistema (13) que contiene r ecuaciones linealmente 
independientes y resuelto con relación а ғ variables. Sin limitar 
la generalidad, se puede considerar que el sistema está resuelto 


* El problema del cálculo de una solución no negativa 
en la que f obtiene su valor máximo зс formula en forma análoga. 
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respecto а las primeras ғ variables de las primeras ғ ecuaciones: 
ха = bı (rr H аваа + + анхо), 
Xa = Dz — (аах, Ж ar+ 2 Xp 42 + — + азах), 


(15) 


En una anotación abreviada el sistema de ecuaciones (15) toma 


la forma 
x= (2 000) i=12. 


Supongamos que los miembros libres b en las ecuaciones (15) son 
no negativos. 

Cada una de las ecuaciones del sistema (15) se puede conceptuar 
como una proyección de la ecuación vectorial 


Ў xP, = Po – (5 2 


en los vectores P,, Ра, ... P,, donde Р, = Р, (1, 0, .., 0); Ра = 
= Р; (0, 1, ... 0) ..; Р, = P, (0, O, ... 1). Los vectores Р, Рз, .., Р, 
forman una base en un espacio r-dimensional. La matriz de 
la descomposición de los vectores Po, Ру, P; ,, m P, en la base 
P., Pa, .., Р, se representa сп la forma siguiente 


(15) 


Py Pz... Р P, Po Pir oe Р, ... Р, 
$ н; bi . а а. 
t2 у Ba ansi - бар . аһ 

: : : 117) 
1 : h amsi . Ay . ûn 
1 b, msi ау. в 


Admitamos que el conjunto de vectores Py, Pa, -~ Р,, Po 1, 
que entran еп la ecuación (16) forma varias bases de un résimo 
orden que se diferencian una de otra aunque sea en un vector, 
Llamaremos a cada una de estas bases base propia. 
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Sean A, da, = А las bases propias de un sistema de vectores 
Pis Pas = Л 5 Preso — Ру, = Pa. Al pasar de una base 
a ctra, todos los coeficientes de las ecuaciones (15) se transformarán 
según las fórmulas (7) del capítulo anterior. Entonces el sistema 
(15) se convertirá en sistemas equivalentes a sí mismo о sea, 
en sistemas con las mismas soluciones. 

Llamaremos la base А positiva con relación al vector cero Po, 
si Po se representa en esta base en la forma 


Po = 3 BIP, Rizo. 


Como se deduce de las ecuaciones (15), la base 4 que contiene 
los vectores P,, Ра, ~, P,, correspondientes a las variables resueltas, 
es positiva, puesto que según la suposición b, > 0. Las variables 
x), respecto а las cuales el sistema (15) está resuelto, 

las restantes, no básicas. La solución de un 
sistema que se obtiene al igualar a cero las variables no básicas 
se llama solución básica. 

En virtud del teorema sobre la descomposición de un vector 
de una sola forma en cada base, la solución básica зе determina 
de una sola forma en cualquier base propia, 

A continuación vamos a estudiar los sistemas de ecuaciones 
algebraicas lineales solamente en bases propias positivas. Por eso 
para reducir la escritura Omitiremos las palabras “propia” 
y “positiva”. 

Examinemos el sistema (15). Definiremos las reglas para pasar 
de una base а otra. Para ello emplearemos la matriz (17). 

1. Encontramos una de las columnas j(r + 1 < j < п) en la cual 
entre los coeficientes аз, азу. ~, ау, ~, а, (1 < i < r) de la varii 
по resuelta ху, hay coeficientes positivos” (en el $ 8 se examina 
el caso en el que no existe una columna tal). 

Supongamos que ésa sea la columna Jy. 


IL Determinamos mín [54 por todos los i para los cuales 
7 
ау, son positivos. 

Supongamos que una de las relaciones mínimas sea aquella 
en la que i = i,. Llamaremos al elemento ai,j, resolvente (o pivote) 
con relación al sistema (15). 


® Se toman en cuenta los signos de los coeficientes 
que se encuentran dentro de los paréntesis en la notación de las 
ecuaciones en forma del sistema (15). 


44 


Ш. Resolvemos la ¡-ésima ecuación con relación a la variable 
х), y la sustituimos por su expresión en las demás ecuaciones del 
sistema (15). Entonces el sistema (15) se transformará en el sistema 
(18) que le es equivalente: 


G=r+ 1, о п), 
y el grupo de vectores Pj, Pz, „Руј, Ру, Pi, 4p, +, Р, se convertirá 
en la nueva base, Ay. Realmente, del (18) se deduce que estos 
vectores forman una base, puesto que la determinante que les 
corresponde es diferente de cero 

Py Pa ... Р, Р, Pre, 


1 


Р, 


o; = ay, >0. 


Esta base es positiva porque los miembros libres nuevos son no 
negativos. Efectivamente, los miembros libres tienen la forma 


bi ь 
Si ay, > 0, entonces según la regla Ш CATA y el miembro 
ма 


libre será по negativo; si ау, < 0, entonces el miembro libre de nuevo 
será no negativo. 
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La transformación del sistema (15) en el sistema (18) determinada 
por las reglas 1, 11, Ш, se llama transformación idéntica O simplex. 

De la definición de transformaciones idénticas se deduce que 
un sistema de ecuaciones algebraicas lineales del tipo (15) puede 
ser sometido cada vez a estas transformaciones si en las partes 
derechas hay coeficientes positivos de las variables. La sucesión 
de јаз transformaciones idénticas de un sistema de ecuaciones 
algebraicas lineales puede ser interpretada Beométricamente como 
una descomposición sucesiva de la ecuación vectorial correspondiente 
а este sistema en las diferentes bases propias positivas que aparecen 
como resultado de una sola sustitución cada vez (о sea en сада 
transformación). 

Teniendo en cuenta que las transformaciones idénticas se emplean 
para la resolución del problema fundamental de la programación 
lineal, las estudiaremos con la siguiente condición complementari 

Supongamos que la última r-ésima ecuación del sistema (15) 
permite una sucesión infinita de transformaciones idénticas de este 
sistema". Eso quiere decir que cualquiera que sea el número 
de transformaciones idénticas que hagamos, en la r-ésima ecuación 
siempre habrá coeficientes positivos de las variables en su parte 
derecha, y en cada transformación, el elemento resolvente no 
pertenecerá a' esta ecuación. 

A la sucesión de transformaciones idénticas le corresponde 
la sucesión А, — Ау = Ay.. de los pasos de una base a otra, 
En el caso de que r y n sean finitos, existe un número finito 
de diferentes bases; entonces en la sucesión Ay => Az — Аз ~.. debe 
haber regreso a bases que ya aparecieron. Diremos que en este 
caso se forma un ciclo cerrado de los pasos de base a base о que 
tiene lugar la “periodicidad”, 

Si en cada transformación idéntica la menor relación de 


b 

min (3, miembros libres b; а los correspondientes coeficientes 
Я 

ар, donde j, es el número de la columna que contiene el elemento 

pivote ay, cs diferente de cero, entonces el valor del miembro 


libre en la r-ésima ecuación del sistema (18) a cada paso sólo 
disminuye, lo que se percibe directamente de la expresión de este 
miembro libre. La elección de la base determina de una sola forma 
el valor de los miembros libres del sistema de ecuaciones, por eso 


| _ "3 Los razonamientos que van a continuación зоп 
justos para cualquier ¿ésima ecuación del sistema (15) (1 <i < 7). 
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a diferentes bases les corresponden diferentes conjuntos de miembros 
libres y viceversa. Por esta razón en el caso dado es imposible 
el regreso a una de las bases que se obtuvieron anteriormente. 
En consecuencia, no puede haber una sucesión infinita de 
transformaciones. Esto quiere decir que después de un número 
finito de transformaciones, о bien todos los coeficientes de las 
incógnitas de la parte derecha de la r-ésimta ecuación se hacen 
no positivos, ay <0 (r + 1 <j < п), o bien el elemento pivote 
en cierta transformación resulta perteneciente a esta ecuación. 


Si comenzando por cierta transformación, el elemento pivote 
perteneciese a ecuaciones con el miembro libre b, nulo, el valor 
del miembro libre en la r-ésima ecuación permaneceria invariable, 
En efecto, en la última ecuación del sistema (18) se ve que en este 


b, 8 
— de la expresión 


ву, 
del miembro libre. En estas condiciones surge la posibilidad де 
la periodicidad. 

Llamaremos hiperplano básico al hiperplano que раза por 
lo menos por dos vectores básicos. 

La existencia de miembros libres nulos en las ecuaciones del 
sistema que se transforma, significa que el vector Ру pertenece 
а uno о a varios hiperplanos básicos. Siguiendo la terminología 
aceptada en la literatura, diremos que en el último caso tiene lugar 
una situación de degeneración. 

Asi que la condición necesaria para que ocurra la periodicidad 
es la existencia de la degeneraci 

Veamos, en calidad de ejemplo, el sistema de ecuaciones: 

xı = 2 — (2x4 + ху), 
х. и 3 — (3 xs) 


Xy = 2 — (ха + 2x5). 


caso es igual а cero cl sustracndo b, — 
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La ecuación vectorial correspondiente tiene la forma xP, + 
+ хаР + xP» = Po — (хаР + xsPs), en la que los componentes 
de los vectores determinan la matriz 
P, P, Р Po Р, Ps 
0.0.2.2 t 
0 1 033 чі 
00121 2 


Los vectores Pi, Ра, Ру, Pa, Ps (fig 17,a) forman varias bases, 
por ejemplo, las bases Pi, Ра, Ру, y Pa, Ps, Pa. А la base Р, 
Pa, Ру le corresponde la solución básica xı = 2; x» = 3; xs = 2; 
xa =0; ха =0. En este caso el vector Ро es una combinación 
lineal con coeficientes positivos 


Ро = 2P, + 3P, + 2Р,. 


Geométricamente, esto significa que el vector Ро “atraviesa” el 
paralelepipedo construido con los vectores P,, Ра, Pa (fig. 17.0). 
Presentemos las mismas ecuaciones en la base Pz, Py, Ра. Para 
ello haremos una transformación idéntica. En la columna de 
la variante х, hay coeficientes positivos. Encontramos el mín 


. Escogemos como elemento pivote al coeficiente de x4 


en la primera ecuación. Resolvemos esta ecuación con relación 
а x4 У Colocamos la expresión obtenida en las otras dos. Entonces 
las ecuaciones toman la forma: 


De las últimas ecuaciones se deduce que en la base Pa, Py, Pa, 
el vector Po ‘se representa en la forma Po = 0P; + Р + Ра. 
En este caso tiene lugar la degeneración ya que el vector Ро se 
encuentra en el plano básico determinado por los vectores Ру y P4 
Por eso Ре está en el lado del paralelepipedo formado por los 
vectores P3, Py, Pa (fig. 17, с). Notemos que si hubiésemos infringido 
las reglas de las transformaciones idénticas, hubiésemos legado 
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a una base en la que Po по podría ser representado como una 
combinación lineal no negativa de vectores de base. Por ejemplo, 


linea, entonces después de las transformaciones, las ecuaciones 
tendrán la forma: 

x= 2— (+2) 

xı = 2 (2x3 — 3х), 

х= —3—(—3ху 7x5). 
De estas ecuaciones se deduce que en la base P,, Pz, Ру, el vector 


Р, se representa como una combinación negativa que contiene 
coeficientes negativos 


Po = -2P, – 3P; + 22, 
y la solución resulta no positiva 


же -2) х= 3; x=0 x=2% x=0 
En este caso el vector Po no atraviesa el paralelepipedo formado 
por los vectores Py, Pa, Pa (fig. 17, d). 
Así que la infracción de las reglas de las transformaciones 
idénticas definidas anteriormente puede llevar a una base negativa, 
Veamos un ejemplo de un sistema de ecuaciones, las transfor- 
maciones idénticas del cual pueden llevar a la periodicidad: 


xı = 0 — (xq = ху = хь + 3x3), 
ма =0 (2x = ки = ха + 


Xy = 1 — (ха + ха + 3x6 — 8x9). 


Examinemos las transformaciones de este sistema con relación 
a la tercera ecuación. En ésta, el coeficiente de х, es positivo 
; 2; L) destacaremos 
ET 

la primera que es una de las menores entre ellas. Por eso se puede 
tomar como elemento pivote el coeficiente de Ja variable x4 de 
la primera ecuación. Después de la resolución de la primera 
ecuación respecto a х, y de la sustitución de x por su expresión 


e igual a uno. Entre las relaciones ( 


49 


en las otras dos ecuaciones obtendremos el sistema transformado: 
ха = Û — (xı — хз + хе + 3x3), 
x=0- (-2 + x + xe - 59) 


(—х + 2xs + 4x6 — Lixa). 


x= 


En la última ecuación en la parte derecha de nuevo hay coeficientes 
positivos de las variables, por ejemplo de xs. El elemento pivote 
en este caso es el coeficiente de xs de la segunda ecuación, puesto 
que el mín = 1 corresponde a la segunda ecuación. Efectuemos 


las transformaciones idénticas que corresponden 
nuevo. Como resultado, obtendremos el sistem: 


х -0-(-2 + х; + sa} 


elemento pivote 


m0 (me + و‎ tp 20) 


ху = 1 — (3x1 — 2ra + ха = хе). 


En la tercera ecuación aparece de nuevo un coeficiente positivo 
en la incógnita хе у se puede tomar como elemento pivote el 
coeficiente de la variable xg en la segunda ecuación. Después 
de la transformación idéntica tendremos: 

xe = 0 — (2x1 + 2x2 + 2х4 — 4x7), 

xs=0-( о ж- а — 3 + ху} 

хум 1—( бху — 4x — 2x4 + 3x7) 


Continuaremos la sucesión de estas transformaciones subrayando 
cada vez en la tercera ecuación uno de los coeficientes positivos 
y destacando con caracteres gruesos su respectivo elemento pivote. 
Obtendremos así el sistema inicial de ecuaciones expresado en 
diferentes bases propias posit 

x =0 — ( xı = 2x — 3x4 + Xs), 

xe = 0 — (2x1 — бх, — 10х„ + 4x5), 

x = 1 (2x, + 2х + Txa — 3x3). 


4-0307 
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РА =0-(-м = 5х, +255 +). 


»=0-( ж+ 2u- xs к), 


ху=1— ( 8x:+17xa—7xs— xe) 


х=0— (хоху x+20) 
x =0 ха xe + 5), 


ху=1- (xu + Xs + 3х, = 8х). 


Según esto, después de seis transformaciones idénticas volvemos 
al sistema inicial de ecuaciones. 

Así, la sucesión de las transformaciones idénticas en las bases 
PPP Pa PPa Ps PaPa PeP, Py PPEP, =P, PPP, Ра Ра 
lleva a la periodicidad. 

No obstante, la existencia de la degeneración no conduce 
obligatoriamente a la periodicidad, lo que se deduce del ejemplo 
siguiente: 


xı =0—( ха—2х + 4x6 + хо = Эха + xo) 
xa 0—(2ха— Xs+ Хе = 2х) — 2ха + 4ху), 
ху = 0—(—х„ + 2x5 — 3х6 + 3x7 — 4х + xo) 
0=1—(3х,— xs 2x + ху Xa+ ху). 


En la última ecuación, el cocficiente de la variable х„ es positivo, 
y como elemento pivote se puede tomar el coeficiente de ха 
en la primera ecuación. 

Efectuando las transformaciones idénticas tendremos: 


ха =0—( хү—2х,+ 4x6 + хә = Эха + ху), 
ха =0—(—2x, + 3x5 — 7х6 — 4х7 + 4хе + 2х0), 
x=0-( х + хе + 4x7 — ха + 2ху), 
0=1—(—Зх, + 5x5 — 14х6 — 2x7 + 8хв — 2х9). 
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Continuando este proceso obtendremos: 
2 1 1 А 4 2 
х=0- (a In E +) 


1 2 2 $ 1 1 
E E E IA з) 


x =0-( х, 7х, + 2х9), 
vtta) 
Y después 
»=0-( т“ + Ја + 4%- Lx + +) 
n=0-( -parpat фос fau- ut фи) 


En la última ecuación el coeficiente positivo de хь es simultánea- 
mente elemento pivote, de modo que no se ha formado una 
sucesión infinita de transformaciones idénticas. Esto muestra que 
la degeneración no ha llevado a la periodicidad. 

Efectuando las transformaciones idénticas que corresponden 
114 


al elemento pivote day = Ê y admitiendo que el valor de las 


variables no básicas es igual a cero, encontraremos una de las 
soluciones básicas del sistema de ecuaciones inicial 


xı =0; x3=0; х= 0; x6=0; 


13 4 1 
ме 96 ie ти 


Es de señalar que la composición de ejemplos con periodicidad 
presenta grandes dificultades puesto que para la realización de 
la periodicidad, los coeficientes del sistema de ecuaciones lineales 
deben satisfacer un gran número de condiciones y solamente una 
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minuciosa selección de los coeficientes puede llevar a la formación 
de un ciclo. 
Mostremos que si se tienen sólo dos ecuaciones 


x1=0- (пао taxt У ao) 
з 
| 09у 
wak (as + axa + У ви) 
Ls 


по puede haber periodicida 


La matriz de los vectores que corresponde 
al sistema (19) tiene la for 


Po Ру Р, Р, 
0 ар ае а) 


1 b an аһ ay 


Si se admite que b> 0, entonces los vectores P, y Р; forman una base, 
Ei valor de los componentes del vector Ру en esta base se puede 
representar con las determinantes 


омеђен ој, а-а 
ау і 0 а 
Sean озу у ај) mayores que сего. Eso quiere decir que se puede someter 
la ecuación (19) a transformaciones, pasando de la base Ру, P а la base 
Pa, Pa y tomando el elemento a, como elemento pivote. En la nueva 
base Fs, Pz, los coeficientes de la ecuación (19) toman el valor determinado 
por las fórmulas (7) del capítulo anterior. Аз, los valores de los 
coeficientes de la variable х, tendrán la forma: 


en la primera línea 


|; 
ы ам і 


en la segunda línea + 


а ا‎ жайды 


вз ам 


аз o 
a if 


Si los numeradores de estas expresiones son positivos, entonces сп la base 
Ру, Pz, los coeficientes de xa serán positivos, Puesto que la relaci 
mínima en el sistema (19) de nuevo corresponde a la primera ecuación, 
tomando el coeficiente de xa de la primera ecuación como elemento 
pivote se puede repetir la transformación del sistema de ecuaciones 
pasando de la base Ру, Pa a la base Pa, Pa. 5 

Supongamos que después de cada transformación se puede realizar 
la siguiente, de tal manera que la sucesión de las bases 


A A A . 
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forma un ciclo. Entonces, cada una de las determinantes, en el orden 
siguiente: д 
(Py, Pa) r Р) (Pei PD _ Фә) (rr) Фара 


St а: a 
(Ру (PP) ... (Pia Pana) (Pam Ps) (Р,Р,) 

debe ser positivo (las determinantes de la fla inferior corresponden 
a la sucesión de los coeficientes positivos en la segunda ecuación, las 
determinantes de la. fila superior, a la sucesión de elementos pivotes) 
En el caso dado los vectores P son bidimensionales y por eso el valor 
de cada una de las determinantes en la sucesión (20) coincide en signo 
y magnitud con el respectivo producto vectorial. Entonces la exigencia 
de que todas las determinantes (20) sean positivas es equivalente a la 
exigencia de que el vector Ру caiga entre los vectores P, y Р; que 
forman un ángulo agudo; el vector Pa, entre los vectores Ру у Ра, etc, 


АЈ 


Continuando. estos razonamientos || mos а que el vector Р, debe 
encontrarse entre los vectores Pa y Р, lo que по puede ser (fig. 18). 
La contradicción también puede ser fácilmente revelada en forma analítica. 
Escribamos el sistema de igualdades vectoriales: 
Ру =091P, + GPa: 
Р, = Ру ЖР, 
Ps = сира аза, 
Pa ¡Pr + ФР, 
Ра аар + аР. 
Moltiplicando vectorialmente cada una de las igualdades (21) una vez 
a la derecha por P, y la otra vez a la izquierda por el segundo vector 
que entra en la parte derecha de cada ecuación, nos convencemos de que 
todos los coeficientes de las igualdades (21) son positivos, puesto que, 
según la suposición, todos los productos vectoriales respectivos son también 
perio». Aif pues, para la ecuación Ру =P + уз, tenemos que 
PsP») = аз (PP) + аза (Р.Р), por consiguiente оз: > 0. Multiplicando 
esta misma ecuación а la izquierda por Р, obtendremos que аз: > 0. 
Después, en la expresión de P, sustituimos а рог su expresión, en la 


en 
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fórmula obtenida sustituimos Ру; por su expresión, etc. Como resultado 
tendremos Р, = үзР, + Тара, Yie > 0; Yu > 0, de donde 


Y 


Р, = 1р, За р, 
Y 


та 


Comparando la igualdad obtenida con la última igualdad del 
sistema (21) 


Р, =ацР, + аР, а> 0; а> 0, 
llegamos a una contradicción. 

Mostraremos que para formar un ciclo, en el caso de cualquier 
número n de ecuaciones, п> 2, hay que efectuar no menos де sel 
iransformaciones idénticas. Puesto que examinamos sólo el caso de una 
sustitución cada vez en las bases, se puede esperar el regreso a la base 

icial después de 2k transformaciones, donde k=1, 2, 3, .. El regreso 
en el caso de k=] es imposible, puesto que las determinantes 
Correspondientes a este caso tienen la forma: 


(PiP; e Pao Ру; (РР, Бро Р); (PaPa ~ В - 2 


(РР, .. В, РД; (P,P; ~ Ру .. PD. 

Las dos últimas determinantes se diferencian en una permutación 
de los vectores P, y Ру por lo que tienen diferentes signos. Eso contradice 
la suposición de su positividad, 

Para demostrar la imposibilidad de la formación de un ciclo cuando 
k = 2 hay que mostrar que ninguno de los caminos de las transformaciones 
idénticas que se presentan más adelante es irrealizable. 


1, (Риба ~ Pj, ~ Ру ~ P) (РАР. P), 


~ P) 


= (P,P, Р)- 
(РР, в) 
= (P,P; Р, 

2 (PiP, Pa = Р) = (P, Pa P)» 
(Р.Р, Р)— 
= (P,P; - P) a 
(P,P: ~ Pa, в); 

3. (АРА — Ph, Р, -- Р) (РР... Р, Р)- 
(РАР, ~ Py, Р)— 


(ФР, ~ Ру ~ Ру Р) 
(PIP, o. Р, ~ Pr, Р). 


El primer camino de transición de base a base es imposible, puesto 
que зе reduce al caso de las dos ecuaciones. a 
segundo camino és posible puesto que las determinantes 
зане о гато tempos de ише pecto que ls dean 
ser las dos positivas. , 
Queda por demostrar la imposibilidad del tercer camino. Para eso 
hay. que convencerse de que la suposición sobre la positividad de todas 
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las determinantes del r-ésimo orden anotadas a continuación lleva a una 
icción: 


sE бєй сз Й r 
' ЫЛ + Џ 
h в 
la а 1 Ч 
* a : 1 ы! 
1 We aja Y 
$ ©. = A E чаа #4 
Ж ом, | р 
x E аһ ; $ t 
Б а, + A 1 
+ .. и 
і а 
Lo 
|! 
i 
tl. 
' 
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ntes tendremos: 


1) ag, > 0: 9 адад, — у, > O; 

2) a, >0: 5) ay, > 0; 

D адаш — баца > 9: y -Ô аа = аа, <0; 
м T) ay, < 0. 


De РАДИ 7) se deduce que ај, < 0. De 2), 5), 6) у 7) se deduce que 
а, < 
Introduciremos las designaciones: 


„= an= -=n 
ay, = ke, а = kin, 
j, а = — 


© k, 9, п, kı, qu son mimeros positivos. 
Entonces de la 3) tenemos 646 — kign? > 0 | qe? – kian? > 0, 
de la 4) tenemos — en +k,ken > O | kk, — 1 > 0, 
de la 6) tenemos qun? — kqe? >0 | yn? — kqe? > 0, 
qe? > кат, 
чат? > kge? 
ааст > Кука дет, 1 > kak. 


Comparando la desigualdad 1 > k,k con la obtenida anteriormente 
как > 1 observamos una contradicción. 

Asi que la periodicidad puede ocurrir solamente después de que 
se efectúen no menos de seis transformaciones. 


$ 8. MÉTODO DE CÓMPUTO DE UNA SOLUCIÓN 
NO NEGATIVA DE UN SISTEMA DE ECUACIONES 
ALGEBRAICAS LINEALES 


Sea dado un sistema de m ecuaciones algebraicas lineales con 
п variables: 


аха + аха + + 01 Xy +. Ф й„х„ = b1, 
a31% + азах + me + وکر + دہ + رکرو‎ = Da, 


ax + арха +. + уху + + арх, = by 
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Sin limitar la generalidad podemos considerar que b, > 0 (esto 
se puede lograr siempre multiplicando las dos partes de las 
ecuaciones respectivas por —1). E 

Presentemos el sistema (13) en la forma 


о-һ-($ 005) i=1,2 т en 


Si en el sistema (22) hay alguna variable que entre sólo en una 
ecuación y el coeficiente de la cual, dentro de los paréntesis, tiene 
el signo «+», entonces se podrá resolver esta ecuación con relación 
a esa variable. 

Supongamos que las ecuaciones del sistema (22) se han resuelto 
con relación a todas las variables tales. Entonces, después de la 
respectiva numeración de las variables, el sistema puede ser escrito 
en la forma: 


X = bi = A A + уху + о. + ац), 


Xy = b, (ааах а Хш»; + + 


O= brei (аин ува + + Фу +. 
‚+ аала) 


о, abreviadamente, 
А 
x = bı — (5 а). 


1=1, 2, m l; 1=1,2, m Yoi 
lo+ro=m; l+k=n b,>0; b,>0 


(23) 


donde 


A cualquier ecuación del sistema (23) que no esté resuelta con 
relación a alguna de las variables la llamaremos O-ecuación. 
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Así que cualquier sistema de ecuaciones algebraicas lineales 
puede ser reducido a la forma (23)”. 

Para encontrar una solución no negativa del sistema (23) lo 
sometercinos a unas transformaciones idénticas que satisfagan las 
condiciones siguientes: 

1. Hallamos la O-ecuación, en la que el miembro libre b, es 
mayor que cero. Si no hay una O-ecuación tal, entonces los valores 
de las variables . 


жеђ; x3=0 ((=1,2 m lo j=lo+1, .., k) 


forman una solución no negativa del sistema (23). Supongamos 
que ésta sea la iésima ecuación. 

2. Señalamos en la i-ésima ecuación el coeficiente positivo ay, +"), 

3. Hallamos el elemento pivote ај, y efectuamos en el sistema 
(23) las transformaciones idénticas correspondientes. 

4. Empleamos la ¡-ésima O-ecuación para las transformaciones 
siguientes del sistema hasta que lo resolvamos о establezcamos que 
el sistema (23) es incompatible * E 
5. Después de resolver la iésima O-ecuación hallamos la 
O-ccuación siguiente en la que el miembro libre sea positivo 
y efectuamos con ella operaciones análogas. 

6. Continuamos este proceso hasta que nos libremos de todas 
las O-ecuaciones. 


Observaciones: 1. El número de ecuaciones en el sistema transformado 
puede resultar menor que т, puesto que al introducir en las demás 
ecuaciones la expresión de la variable resuelta algunas de estas ecuaciones 
se pueden convertir en identidades 0 = 0, las que habrá que excluir del 
tratamiento ulterior. Ў 

If. Supongamos que en el sistema inicial (23) о еп el sistema 
transformado se ha encontrado ecuación resuelta con relación 
а cierta variable, la que tiene el’ miembro libre igual а сего: x = 
= 0 – (ЕВух), y es tal que еп la parte derecha todos sus coeficientes 
de las variables son по negativos. Como el papel de tales variables puede 
ser solamente de cero, para que las x, sean no negativas hay que 
igualarlas a cero en todas las ecuaciones con lo que el sistema examinado 


* Si en el sistema (23) no hay ni una ecuación 
resuelta con relación a alguna variable, entonces la =0 y el sistema (23) 
está compuesto sólo de O-ecuaciones. 

+=) Si en la -ёзіта ecuación, las incógnitas no tienen 
coeficientes positivos, el sistema (23) resulta incompat ya que la 
ecuación 0 = b, — (7 духу, en la que 8, > 0 y todos los ay < 0, no puede 
ser satisfecha con ningún valor no negativo de las variables x, > 0. 
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se simplificará considerablemente. En forma similar excluimos las 
O-ecuaciones que tienen el miembro libre igual а cero y los coeficientes 
diferentes de cero de las variables en la parte derecha de igual signo. 
Tomamos el valor de estas variables en las ecuaciones igual a сего. 


Veamos a qué resultados pueden llevar Jas transformaciones 
idénticas consecutivas que satisfacen las condiciones 1—6. 

a) Puede ocurrir que después de un número finito de trans- 
formaciones idénticas el sistema se libre de las O-ecuaciones. 
Evidentemente, el sistema (23) en сме caso es compatible. 
El conjunto de valores de las variables obtenidas al igualar a las 
variables no básicas a cero y a las básicas a los miembros libres 
en el sistema que no contiene O-ecuaciones, es una solución 
no negativa. 

b) Puede ocurrir que después de un número finito de trans- 
formaciones idénticas se ponga de manifiesto que la O-ecuación 
empleada se convierte en una ecuación de forma 

0 = bi — (Lap 
donde bf > 0, ay < 0 para todas las j. En este caso el sistema 
es incompatible. 

©) Por último, puede ocurrir que el sistema по se libre por 
completo de las O-ecuaciones, y las condiciones del empleo de las 
transformaciones idénticas no se alteren. Como al emplear 
consecutivamente las transformaciones idénticas, el múmero de 
O-ecuaciones по puede aumentar, en este caso, cierta O-ecuación 
poses la propiedad de que ella siempre tendrá por lo menos 
un coeficiente positivo de las variables en la parte derecha pero 
nunca le pertenecerá el elemento pivote. Cambiando esta O-ecuación 


0 = b, – (ба) 


рог la ecuación 


عر 


1 
ع‎ E) 


donde є es un número positivo tan pequeño como se quiera, y, 
estudiando la última ecuación junto con las resueltas veremos que 
en el caso с) ocurre la periodicidad 

Lo dicho anteriormente da la posibilidad de hacer la siguiente 
afirmación: 

Sea cual fuere el sistema de m ecuaciones algebraicas con n 
variables, el empleo consecutivo de transformaciones idénticas permite 
después de un número finito de pasos (excluyendo los casos de 
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periodicidad) determinar si el sistema en la región de los valores 
no negativos de las variables es compatible. En el caso de la 
compatibilidad del sistema, en el último paso de las transformaciones 
se forma cierta solución no negativa. 

Una importante propiedad del método presentado es la 
posibilidad de su empleo para cualquier sistema de ecuaciones 
algebraicas lineales, incluyendo también el sistema linealmente 
dependiente, sin que haga falta separar de antemano las ecuaciones 
linealmente independientes. 

Además, como cualquier sistema de desigualdades lincales al 
añadirle variables no negativas puede ser reducido a un sistema 
de igualdades, este método también es utilizable tanto para 
determinar la compatibilidad de sistemas de desigualdades, como 
para encontrar una de las soluciones no negativas en el caso 
de su compatibilidad. Aquí tampoco es necesario separar al poliedro 
de las soluciones en su “forma pura”. 

En la fig. 19 se representa el esquema de bloques del programa 
para calcular la solución no negativa de un sistema cualquiera 
de ecuaciones algebraicas con ayuda de computadoras electrónicas. 

En la computadora se: introducen los coeficientes de las 
variables en el sistema inicial y los miembros libres. Además, 
al resolver problemas de programación lineal, se introducen los 
coeficientes de la forma lineal que se minimiza. La computadora, 
a base del programa, efectúa la sucesión de transformaciones 
idénticas y comprueba las condiciones de compatibilidad del sistema. 

En el caso de incompatibilidad del sistema, la computadora 
da una señal indicadora y detiene su funcionamiento. Si el sistema 
‘es compatible, la computadora calcula la solución no negativa, 
imprime (si es necesario) y pasa al cálculo de la solución 
óptima. 

Veamos algunos ejemplos. 

Ejemplo. Determinar si es compatible el sistema en la región 
de los valores no negativos de las variables: 


1x, + Зха + З =2, 
1۰x, + 2х; + 4х, LA 
—3x,+3x2+ xa 


0=2—(1-x, + 3х, + 3х3), } 


0 =7 = (1-х, + 2x2 + 4х3), 
0= 


— (3x1 + За + ху). 


el sistema es 
incompatible 


Stop (parada del 
funcionamiento 
de la máquina) 


Elección del primer coeficiente positivo] 
ayr en la primera ión 


Elección de los coeficientes positivos ay» 
len la columna j* entre todas las ecuacio-| 
nes det sistema 


ت 


Cálculo de las relaciones L 
а 
2 


Búsqueda del valor mía Û entre los 
ا‎ 
números —L- 
ap 


Resolución de la ecuación /* con relación! 


a la variable xp 
سے‎ 


Traspaso dela ecuación resuelta al primer 
| "rr del sistema de ecuaciones 


Paso al cálculo de la 
solución óptima 


Salida de la solución 
no negativa 


Registro en la memoria del número j* 
de la variable xp 


Exclusión de la variable ху de todas las 
ecuaciones del sistema 


= 


Verificación de la existencia en el sistema 
Ide identidades tipo 0 =0 y reducción 
[correspondiente del número de ecuaciones 


¡Verificación de la existencia de 0-ecua- 
ciones 


= E == 


м 
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En la primera O-ecuación, el coeficiente positivo de ху, igual 
а uno, es simultáneamente el elemento pivote: 
ху =2- (3x2 + 3х), 
0=5—(–ха + ху, 
0=9 – (12 + 8х). 
En la segunda O-ecuación el coeficiente de ху es positivo; 
el elemento pivote pertenece a la primera ecuación: 


eL (+ -2) 


El sistema es incompatible, puesto que en la última ecuación 
el miembro libre es positivo y los dos coeficientes dentro del 
paréntesis, negativos; por eso, a la segunda ecuación no le satisface 
пі un punto рага el que x, > 0; xa > 0. 

Ejemplo. Hallar una solución no negativa del sistema de 
ecuaciones: 


0 =3 — (2x1 + x + 2x3 + 2x4), 
(1) 0=3-(—2х, + x + 4x3 + 2x4), 
0 = 6 - (2x + 2х +7x + 4x4). 


Efectuamos las transformaciones idénticas del sistema: 


: 
аа Garak) 


(2) 0=6- (2x +6% + 44), 
0 =9 — (3х, + 9x3 + 6x4). 


Una de las soluciones no negativas del sistema es 


a+ ху= 1 х=0; х= 0. 
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El método de cálculo de la solución по negativa de un sistema 
de ecuaciones explicado anteriormente puede servir como uno de los 
procedimientos de definición del rango de una matriz. Efectivamente, sea 
dada una matriz 


ча аз 


Un аз 


ам, фа 4н 


Consideremos sus columnas como los vectores Py, Ру, 


la ecuación 
оњ (ў зт) 
үй 


P. Es evidente que esta сом 
(por ejemplo, xı = 1; хут 0; ји 2, 3, 


la ecuación vectorial o-r-(% ме) en forma de un sistema de т 


P, y escribamos 


en la que haremos Po 
по negat 


n tiene solución 
, n). Representando 


ecuaciones y sometiéndola a transformaciones idénticas nos libraremos 
de las O-ecuaciones, El número de ecuaciones resueltas después de la 
exclusión de las O-ecuaciones determina el rango de la matriz, Como se 
puede demostrar, el número de transformaciones que es necesario efectuar 
también es igual al rango de ја matriz. 

Ejemplo. Hallar el rango de la matriz 


2-4 3 1 
1-2 1-4 

0. 1-1 з 

4 4-4 

Efectuaremos la sucesión de transformaciones identicas: 


بخ 


0= 2 - (Ох, -4х, +3ху+ ха). 

0= 1= (xı = 2a + ху +4x +2ху), 
0=0-( س و‎ ха + م‎ + ху), 
= (Ах, — Ixa + 4x3 — 4х, + 5x9) 


xı = 1= (2x + xg = 4х, + 2x5) 

xy + Ixa = 4х3). 
xa xg + dut xh 
о +Mu-3) 


19 paso 
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xy =0-( 9х, = 4x5) 

xı = 1 — (2% — 13x + 6xsh \ э paso 
0=0-( xy + 12x4 = 3xsh 
0=0-(_ ха + 124 – 3x9) 

ху =0 ~ (12 = 3x3), 

ху =0 = ( 9х, — 4х). 3" paso 
х= 1-х, 

ово 


rad 


Rango de la matri 


$ 9. RESOLUCIÓN DEL PROBLEMA 
DE PROGRAMACIÓN LINEAL 


El problema de programación lineal, como fue dicho anterior- 
mente, puede ser planteado de la siguiente manera: 
se da un sistema de m ecuaciones lineales соп n variables 


аах + 3% + m жару + ахо Da 
ах + axa + + азуу + — + ах, = В, 


a3) 


аах + анха +. + ану + = + аых = bm 
у la forma lineal de esas mismas variables / = cx + Саха +. 
+ суху +. + Cn Xy Hay que hallar entre todas las posibles soluciones 
no negativas 


ба, ха, х) 


del sistema (13) una solución (xî; Xi 
adquiera el menor valor posible. 

+ Para resolver este problema calcularemos una de las soluciones 
no negativas del sistema (13) con el procedimiento explicado 
en el parágrafo anterior. Como resultado, el sistema (13) se 


хр) con la cual f 
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transformará en un sistema equivalente ” 
A + ааа +. + арх +. + ацоф), 
% = bh — (rss ТЕТ 


24 


i=1,2 


т гет r+k<n 


Sustituyendo la expresión de las variables básicas en la forma 
lineal с introduciendo las designaciones respectivas para las 
magnitudes constantes, presentaremos / en forma de una función 
lineal 


s=r-( È e) as 


En concordancia con el método simplex, para hallar la solución 
óptima sometemos el sistema de ecuaciones (24) y la función 
lineal (25) a transformaciones idénticas hasta que se infrinjan las 
condiciones de su realización, o sea, hasta que desaparezcan los 
coeficientes positivos de las variables еп la forma (25). Pero como 
lo expuesto en los dos primeros parágrafos de este capítulo sobre 
la periodicidad se refiere también al problema de programación 
lineal, se puede afirmar lo siguiente: 

El empleo de un número finito de transformaciones idénticas 
permite (excluyendo los casos de periodicidad) encontrar una solución 
no negativa del sistema de ecuaciones con la cual la forma toma 
el valor mínimo. 

Se propuesto varios métodos con el fin de evitar la 
periodicidad. Algunos de estos métodos al efectuar las transforma- 
ciones idénticas, exigen la realización de operaciones complemen- 
tarias cada vez que ocurre la degeneración. Cuando ocurre 
la degeneración, el elemento pivote no se determina en forma 
única. Las operaciones complementarias, еп este caso, en cada 
paso de las transformaciones idénticas, permiten elegir el elemento 
pivote de tal manera que la periodicidad se excluya. El caso 
de la degeneración, como ya se ha indicado, siempre tiene lugar 
cuando el vector P, pertenece a uno o a varios hiperplanos básicos. 


tait 3T Después de una numeración correspondiente de las 
incógnitas. 


5-0307 
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El análisis de la resolución práctica de problemas de 
programación lineal confirma la suposición sobre la pequeña 
probabilidad de la periodicidad, además la composición de ejemplos 
con periodicidad ha encontrado dificultades bastante considerables 
debido a la necesidad de satisfacer un gran número de condiciones 
complementarias. 

Aquellos pocos casos en los que se forma un ciclo pueden ser 
hallados directamente al resolver el problema en forma manual 
(sin computadora). 

Por ejemplo, planteando el problema del cálculo de una solución 
no negativa del sistema 


xı =0 — (ху — Xa — Xs + 3x0), 


mm0- (вал xr +), 


en el que la función lineal 
Ј =3 (x + ха + За — 8x0) 


alcanza su menor valor, y efectuando una sucesión де trans- 
formaciones idénticas como en el ejemplo examinado en pág. 49, 
observamos que hemos llegado al grupo inicial de las variables 
básicas. Si no se presta atención a esto, el proceso de resolución 
puede resultar interminable, al mismo tiempo el valor de la 
función lineal quedará siempre constante. 

Al observar la periodicidad hay que-cambiar el orden de las 
transformaciones idénticas eligiendo otro elemento pivote. Así, en 
el ejemplo examinado; procediendo tal como se explica a continua- 
ción, establecemos que la función lineal f no tiene límite inferior: 


xı = 0 = (ху = ха = ха + 3x6), 


a =0 (ән), 


$ =3 — (ху + х + 3x5 — 8x6), 
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En ei caso dado, el valor de la forma puede ser disminuido 
a cuenta del incremento del valor de x4, pero como en la columna 
de los coeficientes de la variable x, по hay coeficientes positivos, 
о sea, que no hay пі un elemento pivote, esto quiere decir que 
no existe ningún límite para aumentar la variable x. Así la 
función f se hace ilimitada en su parte inferior. 

Al resolver estos problemas en las computadoras electrónicas 
hay que prever en el programa la posibilidad de descubrir el 
regreso a grupos de variables básicas que ya se habían encontrado. 
Sólo en el caso de que se revele la periodicidad, se debe utilizar 
la parte del programa en la que se prevén las opera 
complementarias que excluyen la posibilidad de esta periodi 

La operación complementaria más sencilla es la elección de otro 
elemento pivote en el paso en el que se ha observado la repetición. 

Ejemplo. Hallar el menor valor de la función lineal 


J = Sx, — 10ху + 7x3 — За 


en el conjunto de soluciones no negativas del sistema de ecuaciones; 


= (xı + ху + 7x3 + 2x4), 


Oh — (2x1 = ха + 3ху + 3ш), 


0= 4 — (2х, + 2 + 8ху + х, 


ж =2- (2 + 4% +). 


-(-ж + و35‎ + 5 ха 


23: 

2 

T ха 
КУ аи) 
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» = - (ha), ха = 1 — 0х5), 
%=0 -(-=>). ж =+- (3х3), 


$ = 18 —(-43x5) 


= —18 con х = 0 
жао 
же! 


1 
= 


Ejemplo. Entre las soluciones no negativas del sistema de 
desigualdades: 


—20x, + 12x, — 15x, < 60, 
xı + 2x -3x,<6, 
Зх, + 6х, + 4ху < 12, 
—20x, = 15х + 3x, < 60, 
10x, + 5ха = 2ху > 10, 
бху + Тха + 42x; > 12, 


hallar la solución соп la cual la forma lineal f = x; + х; + ху 
alcanza su valor mínimo. 


Resolución. Introduciendo variables complementarias, conver- 
timos el sistema de desigualdades en el sistema de ecuaciones: 


у = 60 — (—20x, + 12× — 15x3), 
Ya = 6 — (ху + 2х — 3x3), 

уз = 12— (3x; + 6х; + 4х3), 

Ya = 60 — (—20x; — 15х + 3x3), 
0 = 10 — (10x, + 5х; — 2x3 = уз), 
O= 42 — (6x, + 7х + 42x3 — ув) 
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y lo sometemos a transformaciones idénticas: 


193.35 
кык сы и). 


Una de las soluciones no negativas del sistema de desigualdades 
0 


: + Expresemos ahora la 


y fa 


es la solución же. 


forma lineal por medio de las variables secundari: 
1 
AE 

Como se puede advertir en la última expresión de la forma 
lineal, en la región de los valores no negativos de las variables 


1 i 
хан уз – 5 36 a 
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сз ya imposible continuar disminuyéndola. Por eso la primera 
solución ху -+ x =0; ху == ha resultado la óptima y el valor 


correspondiente de la forma lineal es f wa = 2. 


$ 10. SOBRE UN PROBLEMA DE МЇМ-МАХ 


Sean dados п vectores m-dimensionales Ру y el conjunto 
T = () de los números reales ty, j = 1, 2, .., п, que les corresponden. 
Examinemos el conjunto Y = {у} de todas las soluciones posibles 
no negativas y = {ху} de la ecuación 
Р + xP ++ xP, = Po, 


donde Ре es el vector dado en el espacio m-dimensional que no 
es igual a cero. 

Convendremos en designar con x, а las variables x #0 en 
cualquiera de las soluciones y, y a los valores correspondientes 
ty, con fy. 

Con И designaremos el sistema de números 1, que correspon- 
den a las variables ху # 0 en la solución у. 

Acordaremos designar con t, al múmero mayor del sistema де 
números (ty). 

Así tendremos que 


ty}. 
máx (15) 


Presentemos ahora el planteamiento del problema: 
Entre todas las soluciones posibles ye Y hay que hallar una 
solución уда рага la cual el correspondiente 


br = min [mfx (0). 


о bien 
trop = MIR (6) 

Llamaremos óptima y designaremos con yop а la solución y 
(si ella existe) que satisfaga esta condición. Estos problemas se 
Патар problemas de min-max. (Un caso particular de este problema 
es el de transporte por el criterio del tiempo, que se estudia 
detalladamente en el capitulo 1V). 


т 


Si un sistema де т ecuaciones algebraicas lineales con п variables 


У ару= Ы, 1<i<m 1<ј<п, 
A 


que corresponde а la ecuación vectorial У, xP, = Ро, se somete 


a transformaciones idénticas, podremos después de un número 
finito de pasos obtener una solución óptima. 

Para ello es necesario: 

1. Calcular la solución inicial mo negativa del sistema de 
ecuaciones. 

2. Hallar entre las variables básicas x, aquella que corresponde 
al mayor valor de t, igual а tp“. 

Excluir del examen posterior las variables secundarias xy, 
las que t> (рн. 

Encontrar еп la і-ёѕіта línea, dentro del signo Y un 
coeficiente positivo ay y determinar cuál es el elemento pivote 
ау, efectuando después las transformaciones idénticas. 

5: Repitiendo (si eso es necesario) la etapa 4 varias veces, 
sacar la variable x, del conjunto de las variables básicas y exclvirla 
del examen posterior. 

6. Repetir las etapas 2—5 hasta que, en cierto paso, en la 
línea correspondiente del sistema transformado de ecuaciones, todos 
los coeficientes de las variables se hagan no positivos. 

La solución básica obtenida en este paso es la Óptima. 

En efecto, supongamos que después de cierto número de 
transformaciones idénticas, el sistema inicial toma la forma: 


x =b,2 


Supongamos que la variable x, corresponde al mayor valor de 
tp“, Entonces, si b, <0, excluyendo х, del examen posterior, 
obtendremos un sistema incompatible. (Si b, = 0, de acuerdo con 
la observación П de la pág. 58 hay que tomar х, = 0 y también 
ху= 0 para las j para las que а, > 0. Después se continúan 
las transformaciones idénticas). 
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CAPÍTULO Ш 


RESOLUCIÓN DEL PROBLEMA DE TRANSPORTE 
POR EL CRITERIO DEL COSTO 


A continuación se estudia uno de los problemas tipicos de 
la programación lineal, el llamado problema de transporte. 
Al planificar la transportación de mercancías con frecuencia surgen 
cuestiones sobre la organización más racional de los transportes. 
En algunas ocasiones eso significa la realización del cálculo de 
un plan de transportes según el cual los gastos sean los mínimos. 
En otros casos es más importante ganar tiempo y por eso se 
plantea el problema de la elección entre los posibles planes 
de transporte aquel que al realizarlo permita hacer llegar las 
mercancías a su usuario en menos tiempo. 

El primer problema ha recibido la denominación de problema 
de transporte por el criterio del costo: el segundo, problema 
de transporte por el criterio de tiempo. 

El primer problema es un caso particular del problema de 
programación lineal y puede ser resuelto con el método simplex 
el que se ha expuesto en el capítulo anterior. No obstante, las 
particularidades de este problema hacen que se resuclva más 
fácilmente con el método combinatorio. Cuando el número de 
puntos de partida y de destino es pequeño, éste método permite 
resolver semejantes problemas sin emplear computadoras. Cuando 
el número de puntos de partida y de destino ез grande, el 
problema solamente puede ser resuelto con el empleo de las 
computadoras electrónicas. Por ejemplo, un problema de transporte 
con 30 puntos de partida y 40 puntos de destino se resuelve 
en la máquina “Strelá” en 25—30 minutos. 

En el capítulo Ш se expone el método combinatorio de 
resolución del problema de transporte por el criterio de costo 
y se presenta un esquema de bloques simplificado para resolverlo 
en las máquinas computadoras. 


$ 11. PLANTEAMIENTO DEL PROBLEMA 


Se puede formular el problema de transporte por el criterio 
del costo de la siguiente manera: 

Designemos con ay, аз, .., a, la cantidad de unidades de 
mercancías de la misma clase que se encuentran en cada uno 
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de los т puntos de partida, y con bı, bz, .., b, la cantidad 
de unidades de mercancías en cada uno de los п puntos de destino. 
Sea xy la cantidad de unidades de mercancías que se planifica 
transportar desde el ¡ésimo punto de partida hasta el j-ésimo 
punto de destino, y су, el costo del transporte de una unidad 
de mercancias desde el i-ésimo punto de partida hasta el j-ésimo 
punto de destino. 
Sea la cantidad de mercancías que se envía de todos los т 
puntos igual a la cantidad de mercancias que se necesita en los n 
puntos de destino. En tal caso se debe cumplir la condición 


«= Ж 26 


Anotaremos las condiciones de problemas de este tipo en forma 
de tabla 2. Llamaremos solución del problema de transporte 


Tabla 2 


а la matriz X = (xy) de un orden de m x n, cuyos elementos 
no son negativos xy > 0, que satisface las condiciones 


È жута (81,2, m т) об) 


У =b G=1,2 m n) (267 


La condición (26) significa que del i-ésimo punto de partida şe 
saca toda la mércancia; la condición (26°) significa que la necesidad 
del јато punto es satisfecha por completo. El problema se 
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reduce al cálculo de los valores no negativos de xy, con los que 
el costo general de transporte 


C= 2 Lo 
ыы 


sea el menor. 

Empleando la definición del producto escalar de las matrices 
se puede formular de otra manera el problema de transporte por 
el criterio del costo. 

Sea dada la matriz С = (су, donde су son números по 
negativos reales. Hay que hallar entre todas las soluciones Х 
anotadas en forma matriz y que satisfagan las condiciones (26) 
y (26°) una solución con la que el producto escalar (СХ) alcance 
el menor valor. Llamaremos óptima a la solución presentada en 
forma de matriz X que satisface esta condición. 


$ 12. SOLUCIONES BÁSICAS DEL PROBLEMA 
DE TRANSPORTE POR EL CRITERIO DEL COSTO 


Examinemos la matriz dada С = (су) y una matriz-solución 
arbitraria X = (xy) (i = 1, 2, , n). Consideraremos 
que m < n, puesto que en el caso contrario se las podria cambiar 
de Jugar. 

Introduzcamos algunas definiciones. 

Llamaremos celda a un par de números reales ij, conjunto 
a un grupo cualquiera de celdas. A una sucesión de celdas en 
forma iji, Аја, ваја, ву, ~» la llamaremos cadena. Una cadena 
será cerrada si tiene la forma 

idi, didas Ый, =s Ши МА. 
Llamaremos ciclo a cualquier cadena cerrada. Un conjunto se 
llamará cíclico si contiene por lo menos un ciclo, y acíclico si 
no contiene ninguno. 

A cada celda ij le corresponde un y sólo un elemento xy 
de la matrizssolución X, y también un y sólo un elemento су 
de la matriz de costos C. Por eso cualquier conjunto de celdas 
es simultáneamente un conjunto de los respectivos elementos, 
tanto de los ху como de los су. 

'Numeremos los elementos de la cadena cerrada Ө, yendo рог 
ella, por ejemplo, en la dirección de las agujas del reloj. Diremos 
que los elementos de esta cadena que tiene números nones forman 
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la semicadena impar Ө", y los elementos con números pares, 
la semicadena par O”. 
Acordaremos designar la suma de los elementos су de, la 


semicadena impar con ) су, y a la suma de los elementos de la 
са 


semicadena par, con У) cy- 
= 


TEOREMA. Cualquiera que sea la solución X соп un conjunto 
ciclico de elementos ху diferentes de cero, existe una solución Y 
соп un conjunto aciclico de elementos уу tal que (СУ) < (СХ) 
y en la que el número de elementos de la solución diferentes de cero 
es menor que el número de tales elementos de la solución X. 
DEMOSTRACIÓN. Supongamos que los elementos ху diferentes 
de cero de la matriz X forman una cadena cerrada O, . Comparamos 
las sumas У cy у Y су; una de las sumas, por ejemplo Een de- 
a ө 
berá ser tal que Zus Eo 


, Formamos a ak la matriz X = (xy) una nueva matriz 
= (x) cambiando elementos en la cadena O, de la siguiente 
ба, 
Хау = хш, + Xf, ша „= Xy, + ч 
Kin = Xan = XP", wy, =, — 


Aqui iı, iajn =» hh son las celdas de la semicadena impar Ө"; 
баја, бај -» Ша» las celdas de la semicadena par Ө”; хре es el 
elemento menor en la semicadena impar ©". Para todas las demás 
celdas ij tendremos ху = xy. 

Se puede advertir fácilmente que la matriz X’, obtenida a partir 
de la matriz X al transladar рог la cadena O el elemento xpi 
también es una solución del problema de transporte, El número 
de todos los elementos diferentes de cero de la matriz X' se hará 
menor por lo menos en uno. 

Como 


(CX) = (Сх) + E су E cy) ху“, 
СЕ: 


entonces, (CX”) < (СХ). 

Continuando la construcción de tales soluciones después de 
un número finito de pasos llegaremos a la solución Y, fa que по 
tiene cadenas cerradas O de los elementos xy diferentes de cero, 
y (СУ) < (CX). El teorema queda demostrado. 
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COLORARIO. Para encontrar aunque sea sólo una solución óptima 
es suficiente investigar todas las posibles soluciones X con conjuntos 
aciclicos de elementos ху diferentes de cero. 

Cualquier solución X, presentada en forma de matriz, cuyos 
elementos xy diferentes de cero forman un conjunto acíclico, se 
llama básica. 

En general, en una matriz X que represente una solución 
cualquiera, puede haber distinto número de elementos x, diferentes 
de cero. No obstante, independientemente de la solución X que 
tomemos, ella tendrá no menos de n elementos diferentes de cero, 
lo que se deduce de la condición (26°) del $ 11. Por otra parte, 
se cumplirá el teorema siguiente: 

TEOREMA. El número N de clementos xy, diferentes de cero, 
de cualquier solución básica satisface la condición n < N < m+ n — 1, 
en la que т es la cantidad de puntos de partida, y п, la cantidad 
de puntos de destino. 


DEMOSTRACIÓN. Previamente nos convenceremos de que son 
justos dos lemas, Construiremos en un espacio (т + n)-dimensional 
т х п vectores y los pondremos en una concordancia mutua de un 
significado único con el conjunto de todas las celdas ij de la 
matriz С = (cy). Pondremos la celda ij en concordancia con el 
vector Ру que tiene todos los componentes iguales a cero excluyendo 
el i-ésimo y el т + ј-ёзіто, cada uno de los cuales es igual 
a la unidad. Es evidente que a cualquier conjunto de celdas 
le corresponde cierto subconjunto de vectores у viceversa, 
a cualquier subconjunto de vectores le corresponde cierto conjunto 
do celdas. 

LEMA 1. Si el conjunto de vectores (Py) es linealmente depen- 
diente, el conjunto de celdas que le corresponde es cíclico. 

En efecto, supongamos que se sabe que el conjunto de vectores 
(Ру) es linealmente dependiente. Entonces existe una combinación 
lineal de estos vectores que lleva al vector cero q(0, 0, .., 0) y por 
lo menos uno de los coeficientes es diferente de cero. Examinaremos 
sólo los vectores que entran en esta combinación con coeficientes 
diferentes de cero. Sea éste, por ejemplo, el vector P,,,. Entonces 
en la combinación lineal obligatoriamente tiene que entrar por lo 
menos un vector con el mismo índice i, por ejemplo el Рад 
Eso tiene lugar porque en el vector q todos los componentes son 
iguales a cero, pero como en la combinación entra el vector Ру, 
para transformar el ¡-ésimo componente en cero hay que tener 
en la combinación aunque ѕса un vector más con el mismo 
componente diferente de cero. Puesto que ahora entra un vector 


т 


que tiene el componente ja, entonces tiene que entrar por lo menos 
un vector más con el mismo componente, por ejemplo el vector 
P,,  Razonando similarmente construimos la sucesión de los vectores 


Paño Pago Pago Рио = 


Ya que se tiene sólo un número finito de vectores, después 
de cierto múmero tenemos que llegar al vector Py y de él, al vector 
Р. Esta sucesión de vectores corresponde al conjunto de celdas 


id ida ida Ы» <<» ido iji 


que es ciclico. 

LEMA 2. Cualquier conjunto compuesto de т + п celdas ез ciclico. 

Para convencerse de la justeza de este lema es suficiente 
demostrar que cualquiera (т + л) vectores entre los construidos 
son linealmente dependientes, Demostrémoslo. 

Analicemos el vector (т + n)-dimensional. 

m componentes a componemes 
а m +d, +1, 0, +1). 

Este es ortogonal con relación a todos los Py. De aqui se 
deduce que todos los vectores Ру pertenecen al espacio (m + n= 
— 1)-dimensional y por eso cualquiera т + n de ellos son linealmente 
dependientes. 

De la condición (26) del $ 11 y del último lema se deduce que 
n < N < т+ п. En los ejemplos se puede advertir claramente que 
existen unos problemas para los que № = п y otros para los que 
М =m + п — 1. Teniendo esto en cuenta, hacemos la conclusión 
definitiva de que N realmente satisface la condición n < N < m + 
+ n— 1 con lo que se demuestra el teorema, 

De los lemas. que зе acaban de demostrar se deduce que el 
conjunto de vectores {Ру} que corresponde a un conjunto acíclico 
de т+ n = 1 celdas es la base de los vectores Ру construidos 
еп un espacio (т + n)-dimensional. Designaremos соп la letra H 
un conjunto acíclico cualquiera de m + п – 1 celdas. 

Teniendo en cuenta el mutuo y único significado de la relación 
entre todas las posibles bases del espacio (т + n)-dimensional 
y los conjuntos aciclicos de m + n— 1 celdas se pueden expresar 
las propiedades ya conocidas de las bases de los espacios 
multidimensionales aplicadas а los conjuntos acíclicos de т + n — 1 
celdas de la siguiente manera: 

1° propiedad. Sea H, un cierto conjunto acictico de m + n — 1 
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e (JEH, У. Entonces el conjunto H,, obtenido al añadir (j) al 
conjunto Н, contiene un ciclo O y sólo uno. 

2° propiedad. Sca (Г, 7) # (i, j) y (ё, $) € O. Entonces el conjunto 
Ну, que se obtiene de Н, al excluir la celda (Г, j), de nuevo 
es un conjunto асјећсо de т+л—1. 

Los conjuntos H, y Ну sobre los que se trata en las 
propiedades | y 2 se diferencian solamente en una celda. 

Definición. Se llaman conjuntos асісісоѕ de т + n — 1 de una 
sustitución cada vez a dos conjuntos acíclicos de m + n — 1 que 
se diferencian solamente en una celda. 

Examinemos una solución básica arbitraria X = (x,). Según el 
teorema (véase la pág. 75) el número N de elementos ху diferentes 
de cero satisface la condición r < N < m +n — 1. Los demás 
elementos m x n= N de la solución son iguales a cero. 

Construyamos un conjunto acíclico Н de т+п— 1 tal que 
todos los elementos xy diferentes de cero de la solución X se 
encuentren en las celdas de este conjunto. 

Definición. Se llaman ceros seleccionados а los elementos ху 
de la solución X iguales a cero y situados en las celdas del conjunto 
acíclico Н де m+ n~ 1. + 

Definición, Se llama selección al conjunto de elementos no cero 
xy de la solución básica X junto con los ceros seleccionados que 
lo complementan hasta el conjunto acíclico Н de m + n — 1. 

Definición. Se llaman x-seleccionados а los elementos су de la 
matriz C que corresponden a los elementos de la elección X. 

En consecuencia, la selección es la solución básica en la que 
de todos los elementos cero ху = 0 se han escogido aquellos que 
completan el conjunto de elementos no cero хуж 0 hasta un 
conjunto acíclico Н de m+n— 1. Finalmente, llegamos a la 
deducción de que para encontrar aunque sea una solución óptima 
es suficiente investigar todas las posibles selecciones. 


$ 13. ELECCIÓN ÓPTIMA 


solución óptima consiste 


La metodología del cálculo .de 
en que, partiendo de cierta selección inicial, pasamos consecutiva- 
mente a otras selecciones con menor valor del producto escalar 
(CX) y después de un múmero finito de pasos llegamos a la 


* El signo € significa que la celda ij no pertenece 
al conjunto Н. 
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solución óptima. Asi surge la necesidad de formar la primera 
selección. 

Veamos el 
selección: 

Determinamos los elementos de la primera linea de la matriz 
X = (xu). Para ello buscamos y encontramos en la primera linea 
de ја matriz С = (су) el elemento menor. Sea este el elemento су, 
Entonces suponemos que x,y, = mín (a; b). Si a, > bj, encontramos 
en la misma línea el siguiente elemento menor que satisfaga 
а la condición су, > су, Y suponemos Xy, = тица, — ху; by). 
Continuamos estos pasos hasta que no se Satisfaga por completo 
la primera ecuación a, = У, хуу. Si en cierto paso de este proceso 

m 
resulta que el resto de a, es exactamente igual al respectivo b, 
entonces tomando xy, igual а сме resto xy, =b, suponemos 
completamente que el siguiente valor de la variable respectiva 
será хуу, = 0. Después pasamos а la segunda linea, a Ja tercera, 
etc. En las columnas, cuyas ecuaciones están completamente 
satisfechas no se apunta el cero. 

De aquí se deduce que la primera solución, obtenida de esta 
manera, es una selección, puesto que el número de elementos 
marcados es igual а т + п — 1 y entre ellos no hay пі una cadena 
cerrada. 

Veamos otro ejemplo para mayor claridad. Hay que formar 
la primera solución para un problema de transporte determinado 
con іа matriz presentada en la tabla 3. 


utente procedimiento para formar la primera 


Tabla 3 


En la primera línea, el menor de los números (8, 3, 5, 2) es el 
número 2, por eso tomamos хуа = min (10; 15) = 10. Como la 
mercancía del primer punto de partida ya está distribuida por 
completo, pasamos a la distribución de la mercancia del segundo 
punto de partida. Para ello encontramos en la segunda línea 
de la tabla 3 el número menor (4, 1, 6, 7). El número 1 es el menor. 
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Por eso suponemos que x = min (15; 10) = 10. Como en este caso 
queda resto en el segundo punto de partida igual a 15— 10 = 5 
unidades, buscamos en Ja segunda linea el siguiente elemento menor. 
Este será el número 4. Suponemos que xz, = mín (15 — 10; 5) = 5. 
Teniendo en cuenta que las necesidades del primer punto de destino 
están cubiertas por completo y que el resto en el segundo punto 
es igual a 0, encontramos el tercer elemento menor en la segunda 
línea, Tal elemento es el número 6. Por eso suponemos que xz, = 0. 
Ahora pasamos a la distribución de la mercancía del tercer punto 
de partida. Con este fin buscamos entre los números de la tercera 
linea (1, 9, 4, 3) el número menor. Como las necesidades del 
primer punto de destino están satisfechas por completo, no tomamos 
en cuenta en el cálculo el número 1 que se encuentra en el cruce 
de la tercera línea y la primera columna. Continuamos buscando 
el menor entre los números (9, 4, 3). Entre ellos el menor es el 
número 3. Suponemos que ха, = min (25; 15 — 10) = 5. Después 
encontramos en la tercera línea el siguiente número menor que 
es el 4 y suponemos que хуу = 20. 

Es fácil convencerse de que la distribución obtenida es la 
solución. Además la solución formada es selección. Efectivamente, 
el número de elementos clegidos es iguala m + n — 1 = 3 + 4 — 1 =6 
y ellos entre sí no forman ciclos. 

Se puede advertir fácilmente que, por el contrario, cualquier 
elemento no elegido хуу forma un ciclo y sólo uno con los elementos 
de la selección. 

A continuación emplearemos el procedimiento explicado para 
formar las primeras soluciones. 

Supongamos que se tiene la primera selección, la solución Xy 
соп el conjunto H,. Sea el respectivo producto escalar (СХ,) 
igual a Сү. Evaluemos cada elemento су que по entra en Н, por 
medio de la expresión 


Qu= Feun Een 


en la que У cy y У су son las sumas де los elementos de las 
se ۴ 
semicadenas impar y par, respectivamente, de la única cadena 
cerrada que forma cada elemento су con los x-seleccionados (el 
elemento evaluado су se toma сото el primero en la cadena O). 
Señalemos el elemento ciy, al que corresponde la menor evaluación 
Af" y su cadena cerrada con los elementos x-seleccionados. 
Formamos la segunda selección Ху para lo que traspasamos el 


elemento mínimo хре de la primera selección de Ја semicadena 
par a la semicadena impar. En este caso, si en la semicadena par 
hay algunos elementos xy iguales al menor э”, entonces, рага 
mayor certidumbre excluimos del conjunto H, la celda ij con el 
elemento хе que se encuentra la primera en la semicadena раг 
al contornearla en la dirección de las agujas del reloj. En el lugar 
de la celda ij, excluida del conjunto H,, introducimos la celda i,j, 
a la que le corresponde с, У el nuevo elemento de la selección 
Хм, = хре (observemos que а cuenta del desplazamiento del 
elemento xp! por la cadena, en la selección X, pueden encontrarse 
varios elementos сего). Como la selección X se diferencia de la 
selección X, solamente por los cambios en la cadena del elemento 
сыз Entonces es justa la desigualdad 


Eeu- Dey <0, 
© > 


el producto escalar (СХ) será menor que el (СХ,) en el valor 
de 


Са, = Deu] xp". 
«“ © 


Teniendo en cuenta que хра puede resultar igual a cero llegamos 
a la conclusión de que С, > С;. Se puede repetir todo el esquema 
de procedimientos realizado y partiendo de la selección Х; formar 
la selección Xy, etc. Como resultado se formará una sucesión 
de selecciónes Жү, Xa, ..‚ Xs, ~ tal que la respectiva sucesión 
de productos escalares С, > С, 2..2 С, >.. es una función no 
creciente del número de la selección. 

Se llama óptima la selección X, con el conjunto Н, si cada 
uno de los elementos су que no entran en Н, tiene un valor 
negativo A > 0. 

Puesto que para cualquier elemento су que no entra en el 
conjunto H, es justa la desigualdad 


Ley Ecy>o0, 

e S 
el paso, por medio de sustituciones únicas en cada vez, de la 
selección X, a cualquiera otra selección no puede disminuir el valor 


del producto escalar con relación al producto escalar C, corres- 
pondiente а la selección X, *. 


>. А continuación se demostrará que la selección 
óptima es al mismo tiempo la solución óptima. 


6-0307 
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TEOREMA. - Para cualquier matriz С = (су) de elementos reales 
еп la que a, > 0, b, > 0, la sucesión de selecciones Xy, Хз, <> Ху, 
después de un múmero finito de pasos se acaba con una selección 
óptima. 

DEMOSTRACIÓN. 

l-er caso. Si al pasar de una selección а otra, al formar 
la sucesión Xy, Xz, ~ en ninguna de las selecciones aparece 
ху = 0, entonces хре traspasado por la cadena también es diferente 
a cero. En este caso el valor del producto escalar en cada paso 
disminuye. Puesto que cuando m y n son finitos puede haber 
solamente un número finito de diferentes conjuntos Н en сада uno 
de los cuales sc define de una sola forma la selección, entonces 
no puede haber un número de pasos infinito. 

2-0 caso. Si al pasar de una selección a otra, en las nuevas pueden 
aparecer elementos ху con valores cero, surge la incertidumbre 
sobre la calidad de finita de la sucesión Х,, Xz, . Esta incertidumbre 
se basa en que en este caso хре =0 y el valor del producto 
escalar no disminuye. No obstante, esta incertidumbre se elimina 
por completo si simultáneamente al problema de transporte inicial 
se examina un problema con la misma matriz С = (су, pero соп 


ai = aj + E, im ly, 2, m т 
bj = by, 12%. (n— 1), 


donde с es un número positivo suficientemente pequeño. Llamemos 
а este segundo problema el c-problema. Formamos para estos dos 
problemas las primeras selecciones X, y Xi. Entonces como є 
es pequeñísimo, los conjuntos Н, у НЕ coincidirán y en la selección 
Ху no habrá elementos сего xi = 0. Efectivamente, la aparición 
de valores cero de xy en una selección del problema inicial ocurre 
en los casos en los que existen combinaciones de líneas y columnas 
tales que se cumple la igualdad 


Хан Yo 
Por eso la aparición de valores cero ху en el problema puede 
suceder en dos casos: 


а) у а +а6= Ўы, 


iah Jah 
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29 > 
b) Ба+ае= Y 5) + b, + те. 
iak ish 


Escogiendo e diferente de las raices de las ecuaciones expresadas 
en las condiciones a) y b) excluimos las posibilidades de que 
aparezcan elementos cero en la selección Xi. Dado que la matriz 
С = (су) es la misma para los dos problemas, coinciden los 
elementos Су con la mayor evaluación negativa, las cadenas cerradas 
de estos elementos con los elementos de las selecciones X, y Xf, 
así como la situación en los conjuntos Н, у Hý de los elementos 
x" y хўр, destinados a ser transladados por la cadena. El paso 
de la selección X, а la selección Ху en el problema inicial es 
de hecho el paso de la selección Xî а la selección X$ en el 
e-problema. Éstos razonamientos, teniendo en cuenta que є es 
pequeñísimo, son correctos en cualquier paso. Como el e-problema 
es un problema en el que se carece de xy = 0 cuando tiene lugar 
la disminución del producto escalar en cada paso, pues en este 
s-problema después de un número finito de pasos llegaremos 
a la selección óptima Xf. Debido a la coincidencia de las matrices 
С = (cu) y los conjuntos H, y Hi, la selección X, también es 
óptima. 
El teorema queda demostrado. 


$ 14, INVARIABILIDAD DE LA SUCESIÓN 
DE SELECCIONES, EQUIVALENTES A LAS TRANSFORMACIONES 
DE LA MATRIZ DE LOS COSTOS 


El cálculo de los elementos está ligado a un gran volumen 
de cómputo. Surge la tarea de definir transformaciones suficiente- 
mente sencillas de la matriz inicial С = (с), que sin cambiar 
la sucesión de selecciones, reduzcan considerablemente el proceso 
de búsqueda de la selección óptima. 

Para ello introduzcamos el concepto de transformación equi- 
valente de una matriz. 

Definición. Sea dada la matriz С = (су) y los números arbitrarios 
тїт е} Fmi эы S23 mi зе Diremos que la matriz D=(dy) 
es equivalente а la matriz С = (су), si ella ha sido obtenida de 
la matriz С mediante la fórmula D = (су + r; + sj). 

TEOREMA. La sucesión de selecciones es invariante con relación 
a las transformaciones equivalentes de la matriz. 
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DEMOSTRACIÓN. Transformemos la matriz C en una triz D 
equivalente. Formamos para la matriz C la selección inicial Xy 
y partiendo de ella, por medio de una sola sustitución cada vez, 
formamos la sucesión Ху, Xz, Ху, .., que converge hacia. la 
sciección óptima Ху. 

Partiendo de la misma selección X,, formamos para la matriz 
D la sucesión Ху, Ха, X3, .., convergente hacia la selección 
óptima Xi. Fácilmente se puede ver que la diferencia de las 
evaluaciones Дб — Ag para cualesquiera dos elementos сб y ср 
que no entren en la selección X, es igual a cero: 


Af = Ab = [Z e “E - [Z(G ++ s) 
-Eleh +n + s)) = [E9 – E – [E (e) – 260) – 
= ES 7 


EM +s) + (+) =0 
= е 


Las últimas dos sumas son iguales entre sí y tienen diferente signo 
puesto que cualquier número ғ (ó з) contenido en la semicadena 
se contiene también en la semicadena par. Como la 
¡ón de los elementos en las matrices С у D coincide, 
coinciden también los elementos que tienen el menor valor. Por 
eso, al pasar de la selección X, a la selección siguiente tanto 
para la matriz С como para la matriz D, las selecciones X, у Ху 
coinciden, Estos razonamientos son justos en cualquier paso lo que 
quiere decir que la sucesión de selecciones Xi, Xa, .., Xy no 
cambiará al realizar transformaciones idénticas en la matriz. 
El teorema queda demostrado. 

Es evidente que esta afirmación no se altera al repetir las 
transformaciones idénticas en la matriz C. 

Teniendo en cuenta que los elementos de cualquier selección 
no forman cadenas cerradas entre sí, en cada paso se puede 
realizar una transformación idéntica que reduce a cero los elementos 
x-seleccionados. Para ello, por ejemplo, es suficiente añadir а los 
elementos de cada columna (línea) un número con el signo 
contrario al de suma algebraica del número que se añade a la 
línea (columna) y el valor del elemento x-seleccionado que se 
reduce a cero. 

La matriz con elementos cero x-seleccionados tiene la propiedad 
de que la evaluación de cualquiera de sus elementos que no sean 
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los x-seleccionados será igual al mismo elemento. Por eso по cs 
necesario hacer grandes cálculos para determinar las evaluaciones 
de todos los elementos. Para formar la selección siguiente es 
suficiente marcar el elemento negativo mayor (según su valor 
absoluto). Si ya todos los clementos de la matriz transformada 
resultan no negativos al mismo tiempo que los elementos 
x-seleccionados son cero, entonces la última selección es ya la 
óptima. De la invariabilidad de la sucesión de las matrices se 
deduce que la selección óptima es al mismo tiempo la solución 
óptima. Efectivamente, la solución óptima con la matriz trans- 
formada de manera equivalente, en la que los elementos 
x-seleccionados son iguales a cero, tiene la propiedad de que su 
producto escalar es menor (o por lo menos igual) que el producto 
escalar de cualquiera de las soluciones. Por csa razón la selección 
óptima que hemos encontrado es la solución óptima. 


$ 15. ALGORITMO DEL CÁLCULO DE LA SOLUCIÓN ÓPTIMA 


J. Los datos del problema se anotan en tablas. 

П. Se define la primera selección. 

TIL. Se reducen los elementos x-seleccionados a cero, Si después 
de eso los demás elementos de la tabla transformada resultan 
no negativos, entonces la primera selección ya es la óptima. 

IV. Si después de la reducción а cero de los elementos 
x-seleccionados en la tabla hay números negativos, encontramos 
el elemento negativo de mayor valor absoluto. 

У. Formamos una cadena cerrada del elemento negativo mayor 
con los elementos x-seleccionados reducidos a cero y componemos 
la segunda selección transpasando por la cadena el número Xy, 
el que sea el menor en la semicadena par. 

VI. Сото el elemento negativo mayor ha resultado un 
x-seleccionado, lo reducimos a сего, pero de tal manera que queden 
iguales а cero los demás elementos x-seleccionados que corres- 
ponden a la nueva selección. 

УП. Repetimos estos pasos hasta que lleguemos а una selección 
tal en la que en la tabla transformada todos los elementos 
x-seleccionados sean iguales a cero y los demás sean no negativos. 
Esta última selección será la solución óptima. 

УШ. Para calcular los gastos de transporte, que corresponde 
a la solución óptima, se deben reunir en una misma tabla la matriz 
inicial de costos y la última solución. La suma de los productos 
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de los dos números que se encuentran en cada una de las celdas 
de la tabla determina el valor de los gastos de transporte. 

Para mayor daridad, se presenta la interpretación geométrica 
del proceso de cálculo de la solución óptima (fig. 20). 


Ст) - costo iniciat 
са) сова final 
o 


Marquemos en el eje horizontal los números de las soluciones 
y en el eje vertical, los respectivos valores de. los gastos. 

Al comenzar, formamos la primera solución y encontramos uno 
de los puntos, el Ај. La reducción а cero de los elementos 
x-seleccionados nos Пеуа-ај punto 4, del eje horizontal. Si en 
este punto la tabla transformada no tiene elementos negativos, 
la solución que corresponde al punto A, es la óptima. Si la tabla 
transformada no tiene esta propiedad, encontramos el punto A3; 
de nuevo volvemos al eje en el punto А; y comprobamos si hemos 
alcanzado уа la solución óptima. Si по la hemos logrado todaví: 
encontramos el punto Ay, etc. 

Después de un número finito de pasos se llegará forzosamente 
al último punto, el que corresponde a la solución óptima. 
Obtendremos la linea quebrada de los gastos que aparece en la 
fig. 20. Como se puede ver en el dibujo, después de cada paso 
el valor del gasto disminuye (o en todo caso no aumenta). 

Continuemos la resolución del problema de transporte planteado 
en la pág. 79. En la tabla 4 se presenta la primera selección para 
este problema. Reduzcamos а cero los elementos x-selecciónados. 
Para ello restemos 4 de los elementos de la primera columna 
de la matriz de costos; de los elementos de la segunda columna 
restemos 1. Como resultado de esto se reducirán a cero los 
elementos xrseleccionados que se encuentran en los cruces de la 
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segunda línea con la primera y la segunda columna, Restemos 6 
de los elementos de la tercera columna, entonces el elemento 
s-scleccionado que se encuentra en el cruce de la segunda línea 
соп la tercera columna también se reducirá a cero y el elemento 
de la matriz de costos que se encuentra en intersección de 
la tercera línea con la tercera columna tomará un valor igual 
а (4— 6 = —2). Añadamos а los elementos de la tercera línea 
de la matriz de costos el número 2 у de los elementos de la 
cuarta columna restemos el número 5. Entonces, los elementos que 
Se encuentran en los cruces de la tercera línea con la tercera 
y la cuarta columnas se reducirán a cero. 


Tabla 4 


Л УТ 1 
ОАО О 7 | 
ZA TADA 


Por último, queda reducir a cero el elemento x-seleccionado 
que se encuentra en el cruce de la primera línea con la cuarta 
columna. Para ello hay que añadir a los elementos de la primera 
línea el número 3. Como resultado de esta transformación idéntica 
de la matriz de costos, la tabla 4 se convertirá en la tabla 5. 

Como se ve en la tabla 5, hay un elemento negativo (—1) 
Por eso no podemos estar seguros de que la primera selección 
es la óptima. Como en este ejemplo hay un solo elemento negativo, 
el mismo es el de mayor valor absoluto. 

Formemos una cadena cerrada de este elemento con los 
x-seleccionados у compongamos la segunda selección. Eso se 
muestra en la tabla 6. Con este fin transpasemos de la semicadena 
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par а la impar el número de unidades de mercancías que en la 
semicadena par sea el menor. 

Como el número —1 уа es xeseleccionado, hay que reducirlo 
a cero, para lo cual sometemos de nuevo la matriz de costos 
a transformaciones idénticas. Añadimos a los elementos de la 
primera columna de la tabla 6 un múmero igual a 1. Obtendremos 
como resultado la tabla 7 en la que no hay ni un elemento 


negativo, y todos los elementos x-seleccionados son iguales а cero. 
Tabla 6 


La selección presentada en la tabla 7 es, basándonos en lo 
indicado, la solución óptima. Reunimos en una sola tabla 8 
la matriz de costos inicial y la solución óptima. El gasto de 
transporte que corresponde а este plan es C=2-10+5-6+ 
+ 101 + 5<3 + 4-15 + 1-5 = 140 unidades de costo. 

Según este plan hay que organizar el transporte de la 
siguiente manera (fig. 21): Toda la mercancía del primer punto 
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de partida se transporta al cuarto consumidor; el segundo punto. 
de partida abastece con 10 unidades de mercancia al segundo 
consumidor y con 5, al tercero; desde el tercer punto de partida 
se tansportan 5 unidades de mercancias al primer consumidor, 
15 al tercero y las 5 unidades de mercancías restantes se envían 
al cuarto consumidor. Con este plan, el gasto de transporte 
(a igualdad de las demás condiciones) será el mínimo. 


Examinemos un ejemplo más complicado. 

Ejemplo. Hallar el plan de transporte óptimo de una mercancía 
homogénea desde cuatro puntos de partida hasta seis puntos 
de destino para las condiciones definidas en la tabla 9. En esta 
tabla a, у фу se miden en miles de toneladas y су, en miles 
de rublos. 
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Escribamos otra vez la tabla inicial, introduzcamos en la parte 
inferior de las celdas los números de la selección inicial y presentemos 
todo el proceso de obtención de la solución óptima como una 
sucesión de las tablas 9а, b, с, d, e, f. 

Tabla 9а 


Primera 
selección 


C,=75 Е 
de la lablo 


Traspasamos 
$ unidades por 
la Cadena. 

El gasto dis= 
minuyo en 

20 unidades 


Tabla 9с 


El elemento que se oscluye de la 
Solución persed" sombresdo 
Segunda Transposamos 
selección 0 unidades por 
Ca 55 la. cadena. 
El gasto queda 
el mismo 


Tronspasamos рог 


AAAA VAA и 


оь l A a RE 
6255 AAA en 2х2+2х1=6 


ATA VA AN 
ARA 


e 


ARNS 


4 
ЛУТ жинге 


unidad рог cadeno. 
El gosto disminuye 
en Ixi=1 unidad 


Unamos la tabla inicia! de costos 
Y la solución óptima (tabla 10). 


En la última tabla se indica cómo hay que organizar el 
transporte para garantizar el mínimo de gastos. La variante 
óptima en comparación con la primera variante garantiza una 
disminución de los gastos de transporte de 27 mil rublos, ya que 
estos gastos en la primera variante son de 75 mil rublos y en 
la óptima sólo de 48 mil rublos. 


Es de interés observar que la tabla que contiene la variante 
óptima (la 5* selección) se indica по sólo la única solución 
óptima, sino todas las soluciones con los mismos gastos. 

Efectivamente, en esta tabla hay elementos cero que no están 
marcados con un círculo, Marquémoslos con una estrella. Cada 
uno de estos elementos forma una cadena cerrada con los elementos 
de la solución. 

Empleando con las reglas anteriores el procedimiento de una 
sola sustitución cada vez se pueden obtener otras soluciones que 
tienen los mismos gastos. 

Esto permite tomar como óptima una solución, por ejemplo, 
como la presentada en la tabla 11. Los gastos correspondientes 
а este plan de transporte también son iguales a С=1.2+4.3+ 
+4:1+1.5+1-1+1-2+1-2+5-.1+3-5 = 48. 

Fácilmente se ve que en este caso cada cadena cerrada compuesta 
sólo de elementos cero, tiene la propiedad de que son iguales las 
sumas de los gastos en las semicadenas par e impar de la 
respectiva cedena cerrada de la tabla inicial lo que permite 
traspasar xpi por la cadena sin alterar el gasto total. En una 
serie de casos resulta muy útil definir todas las soluciones óptimas. 
Como veremos más adelante, aplicando al conjunto de soluciones 
óptimas el procedimiento de elección de la variante óptima de 


Tabla 11 
у> 
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transporte bajo el criterio del tiempo, se puede encontrar un plan 
de transporte que se realice en un tiempo mínimo. 

Por último, refirámonos a la composición de la cadena cerrada 
formada por el elemento negativo mayor con los elementos 
3x-seleccionados. Cuando m y п no son muy grandes, las cadenas 
se ven directamente como en los ejemplos examinados anterior- 
mente. En cambio, si m y n son relativamente grandes, la búsqueda 
de la cadena cerrada se realiza de la forma más sencilla de la 
manera siguiente: se tachan las líneas y las columnas que contienen 
un cero marcado (puesto que para componer la cadena hay que 
tener no menos de dos ceros en la línea o en la columna); 
después de tachar la primera, se tacha la segunda, etc., hasta que 
lleguemos a una situación en ја que enseguida se vea la cadena 
cerrada. 

Por cuanto la resolución detallada de un ejemplo que tenga 
una tabla de grandes dimensiones nos ocuparía mucho lugar, aquí 
nos limitaremos a examinar un caso еп el que se debe hall la 
cadena cerrada. Esto se presenta en la tabla 12. Al principio 
tachamos todas las columnas que contienen un cero (claro está, 
excluyendo la columna en la que se encuentra el elemento negativo 
mayor). Después examinemos todas las líneas y tachamos aquellas 
en las que ha quedado un cero. De nuevo revisamos las columnas 
y hacemos lo mismo. 

Ahora para componer la cadena cerrada solamente hay que 
unir los elementos que quedan comenzando por el número negativo 


Tabla 12 
алеу ажа aa 
РТ x т), 
LA я 2 
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y revisando los elementos cero, por ejemplo, еп Ла dirección 
de las agujas del reloj. 

Hemos estudiado la metodología de resolución de los problemas 
de transporte Бајо el criterio del costo a condición de que 


У а= Y Ы. En este caso, la cantidad de mercancia en todos 
ата 
los puntos de partida es igual а la cantidad de mercancia necesaria 
en todos los puntos de destino. 

Si 


2a> уь, 
4 کر‎ 


Tabla 13 Tabla 14 


entonces introducimos un punto de destino ficticio con un consumo 
һа=ўа- У 
MM A 


y соп шп costo de transporte en este punto €, „=, =0 para 
їе 1, 2, ‚т 

Ejemplo. En cuatro depósitos hay respectivamente 70, 40, 50 
y 20 toneladas de combustible. Se debe planificar el transporte 
de combustible a tres consumidores de tal manera que el gasto 
del transporte sea mínimo. En la tabla 13 se dan las condiciones 
del problema. En este caso es preferible formar la primera 
solución por las coluńnas. Esta solución está marcada соп 
círculos en la tabla 14. 

Es fácil comprobar que la solución presentada en la tabla 14 
es la óptima. De acuerdo con esta solución al primer consumidor 
hay que llevarle 10 toneladas del primer depósito y 20 toneladas 
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del cuarto; al segundo consumidor, 30 toneladas del primer 
depósito y 10 toneladas del segundo; al tercer consumidor, 
30 toneladas del segundo depósito; además en el primer depósito 
quedan aún 30 toneladas de combustible; en estas condiciones 
no es conveniente emplear las reservas del tercer depósito. 


Composición de la selección inicial 


Reducción a cero de los elementos x-se- 
leccionados de la matriz de los costos 
(transformación de la matriz) 


Búsqueda en la matriz transformada del 
elemento negativo de mayor valor absoluto! 


Control del final de la resolución: se com- 
prueba la existencia de elementos negati- 
vos en la matriz transformada 


¡Cálculo del costo del 
transporte según el plan 
óptimo 


Salida a la impresora 
[del plan de transporte 


[Composición de la cadena cerrada formada | . 
óptimo y el valor del 


рог el elemento negativo de mayor valor 
absoluto y los elementos x-seleccionados | | costo correspondiente 


1 1 


Composición de una nueva sección | | Stop 


AA 


mn 


El algoritmo expuesto para la resolución del problema de 
transporte bajo el criterio del costo puede presentarse como un 
proceso rutinario de operaciones aritméticas y lógicas. Este proceso 
se realiza fácilmente en una máquina computadora electrónica. 
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En la fig. 22 se ilustra un esquema de bloques para la resolución 
de este problema. A Базе de este esquema de bloques fue preparado 
un programa para la resolución del problema en la computadora 
“Strel: 

Para una matriz de un orden de m x и < 500 el tiempo de 
resolución del problema no es mayor de 8—10 minutos; para una 
matriz de un orden de m x n < 1500 pueden hacer falta para resolver 
un problema hasta 25—40 minutos; para una matriz de un orden 
mx n de varios miles hacen falta hasta 8—10 horas para resolver 
un problema en la computadora “Strela”. 

Por ejemplo, en la computadora “Strelá” se resolvió el siguiente 
problema de transporte. 

Desde nueve puntos de partida hay que transportar una 
mercancía a catorce puntos de destino (tabla 15). La selección 
inicial obtenida en la computadora está representada en la tabla 16. 


Tabla 15 
MAREAS " 0 
(Yue Yo Te (У To Te To Too] 


els Гау Терге] 


' 2 22 567 • ЖА ones 


DEA: 


La selección final que corresponde а la solución óptima se ilustra 
en la tabla 17. 

Es fácil calcular que los gastos que corresponden a la selección 
inicial son de un valor igual а С = 971 у los que corresponden 
a la óptima, С = 703. 

El tiempo de resolución de este ejemplo (en un solo cálculo) 
fue de menos de un minuto. 

Otro problema de transporte con una matriz de costos de 
т x n = 30 х 38 (solucionada en dos cálculos) fue resuelto en un 
lapso de 37 minutos. 


7.0307 
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CAPÍTULO IV 


RESOLUCIÓN DEL PROBLEMA DE TRANSPORTE 
POR EL CRITERIO DEL TIEMPO 


La economia de tiempo, al elaborar un plan de transporte, es 
absolutamente necesaria en una serie de casos prácticos importantes. 
Por ejemplo, al transportar productos de fácil deterioro es necesario 
su traslado a los puntos de destino en el plazo mínimo posible, 
En la época de la cosecha un factor muy importante es el rápido 
traslado del grano a los centros de acopiamiento. 

Los problemas de'este tipo son los de transporte por el criterio 
del tiempo. A continuación se estudia uno de los algoritmos de 
la resolución del problema de transporte por el criterio del tiempo. 


5 16. PLANTEAMIENTO Y RESOLUCIÓN DEL PROBLEMA 


Sean dados m puntos de partida de una mercancía homogénea 
y п puntos de destino. Sean as, аз, m а, ..‚ ат la cantidad 
de unidades de mercancía que se encuentran €n al primero, segun- 
do, ... Fésimo, .., mrésimo puntos de partida y bi, bz, « бр => Da 
la cantidad de unidades de mercancia que tienen que ser entregadas 
en cada uno de los n puntos de destino. Sea ty el tiempo 
(en días o en horas) necesario para el transporte de la mercancía 
desde el i-ésimo punto de partida hasta el j-ésimo punto de destino. 
Sea xy la cantidad de unidades de mercancía que planificamos 
transladar del iésimo punto de partida al j-ésimo punto de 
destino. Hay que hallar el plan de transporte óptimo, es decir, 
unos números no negativos xy don los que el tiempo de transporte 
de todas las mercancías necesarias a todos los puntos de destino 
sea el mínimo. 

El planteamiento matemático de este problema se reduce 
a lo siguient 

Se da un sistema de ecuaciones algebraicas lineales: 


A xy =a, ((=1,2.,m, 
Ф qn 
РЕ 0=1, 2... л), 
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en las que los números a, y b; satisfacen la condición 
Na Ed 
дч Pr 


Sea, además, dada la matriz del tiempo Т = (tų). A cada solución 
по negativa X = (хц) del sistema de ecuaciones (27), es decir, 
a cada plan de transporte le corresponde una determinada matriz 
Ту = (tf), cuyos elementos son (f = ty si xy > 0 у) =0 si xy =0. 
Llamemos пр“ el elemento mayor de la matriz Ту. 

Se debe hallar una solución no negativa Хад para la cual 
ври, Sea el menor entre todos los гр" que corresponden а diferentes 
soluciones no negativas X. Es cómodo presentar las condiciones 
de сме problema en forma de la tabla 18. 


Tabla 18 


Como se puede mostrar fácilmente, este problema no se resuelve 
con la ayuda del algoritmo que fue descrito para el problema 
de transporte por el criterio de gastos. 

Por ejemplo, la solución para la que la forma lineal С = J tyxı, 
alcanza el valor mínimo, como regla general, no es la óptima por 
el tiempo. Para el ejemplo dado en la tabla 19 e ilustrado en 
la fig. 21 se presenta la solución óptima por el mínimo de la 
forma lineal С = Ууу en la tabla 20. En la tabla se ve que toda 
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la mercancia será trasladada hasta los puntos de destino solamente 
después de seis días. Por otra parte, la variante de transportes 
presentada en la tabla 21 corresponde a un cierto valor mayor 
de la forma lincal, pero el traslado de las mercancías a los puntos 
de destino se cumple en cuatro días. 

Al transportar productos de fácil deterioro los gastos en el 
recorrido adicional de vagones (camiones) se cubren con la 
conservación de la. calidad de miles de toneladas de productos 
alimenticios destinados para la población. 

La demostración del método de resolución que se describe más 
adelante se deduce de la teoría general de la programación lineal 
que fue expuesta en'el capítulo H. 

Veamos la descripción general del método en la resolución del 
problema que se plantea con la tabla 22. En esta tabla se ve que 
hay que transportar 125 unidades de mercancías que se encuentran 
en seis puntos de partida a siete puntos de destino. Supongamos 
que los números de esta tabla significan eb tiempo necesario, 
en horas, para el transporte de la carga desde el respectivo punto 
de partida hasta los de destino. Por ejemplo, el número 31 que 
se encuentra en el cruce de la cuarta linea con la cuarta columna 


Tabla 21 


Tabla 22 
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muestra que el tiempo de transporte de la carga desde el cuarto 
punto de partida hasta el cuarto punto de destino es 31 horas. 
Hay que hallar el plan de transporte óptimo por el tiempo. 

Representemos la tabla anotada (llamémosla la tabla pequeña) 
en forma de una tabla grande (tabla 23). 

La línea superior contiene los índices de las variables x; desde 11 
hasta 67, la inferior, los valores de ty contenidos en las respectivas 
celdas ij. Por ejemplo, el par de números 43 que están en la línea 
de arriba significa el lugar de la variable x.y, el valor del tiempo 
tas que corresponde a esta variable es 39; рог eso 43 y 39 se 
encuentran en una columna. 

Toda la tabla está dividida en zonas. El número de zonas es 
igual al número de puntos de partida o, lo que es lo mismo, 
al número de líneas en la tabla inicial. Además, la tabla está 
jvidida horizontalmente. En la mitad superior de las zonas hay 
líneas: el número de líneas en la tabla inicial; en la parte 
inferior, hay siete líneas: el número de sus columnas. En la 
semizona superior al nivel de cada número figuran siete signos 
menos en cada respectiva zona vertical. Por ejemplo, al nivel del 
30, los menos están anotados en la cuarta zona. Si examinamos 
la tabla inicial veremos que esta anotación significa la expresión 
del resto de la carga en el cuarto punto de parti 


30— хар — Хаз — Хаз — X44 — Xas — Xas — Хал. 


Similarmente a esto, en la parte inferior la colocación de los signos 
menos en cada línea expresa la falta de mercancía en los puntos 
de destino. Por ejemplo, la colocación de los menos en la línea 
donde está el 27 indica que la falta de mercancía en el punto 
de destino que corresponde a la cuarta columna es igual a 


27 — Xia = Ха — Хуа — Kaa ¬ Xg4 — Хз. 


Con otras palabras, la colocación de los menos en la tabla 
corresponde a las ecuaciones que expresan las condiciones del 
problema. 

Construyamos la primera solución utilizando el procedimiento 
empleado en el capítulo Ш. Distribuyamos las 15 unidades del 
primer punto de partida entre los puntos de destino hasta los que 
el tiempo del recorrido es el menor; las 7 unidades del segundo 
punto de partida, entre los puntos de destino que todavía no 
están abastecidos hasta los que el tiempo del recorrido es el menor, 
etc, etc. Así se forma la primera solución presentada ел la 
tabla 24 que, en general, puede no ser la óptima. Como se ve, 
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realización exige 28 horas para el transporte de la mercancía. 
Eso nos permite, a continuación, excluir de la revisión las celdas 
que tienen ty > 28, lo que se ha hecho en la tabla anexa. 

Para la composición de la siguiente variante de solución, cuya 
realización exija menos де 28 horas, dirijámonos a la tabla grande, 
Para los datos de esta tabla, la composición de la selección se 
hace de la siguiente manera: en la línea inferior de los valores ty 
se halla el menor de los números de la primera zona. Ese es el 
número 7 que corresponde а la columna 14, o sea, a la variable x4. 
Examinando esta columna nos convencemos de que podemos hacer 
ха igual а 15, ya que la reserva de mercancia es de 15 unidades 
y la necesidad, de 27. Eso permite anotar la variable x,4 а la 
izquierda de la igualdad en la primera línea y suprimir el signo 
menos que se encuentra en el cruce de la primera línea y la 
columna 14, El segundo menos de esta columna que se encuentra 
al nivel del número 27, también se puede suprimir, y de la línea 
que corresponde al número 27 se resta la línea que corresponde 
al número 15. 

En las líneas primera y décima de la primera zona han habido 
cambios. Estas líneas tienen ahora la siguiente forma: 


ха = 151—1—1—11—1— 


Pasamos a la segunda línea y a la segunda zona. Hallamos 
el menor de los números ty de la segunda zona. Ese es el número 
6 que corresponde a la variable ха. Revisando la columna 
verticalmente, elegimos el menor de los доз números que se 
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encuentran a la izquierda al nivel de los signos menos de esta 
columna. Ese es el número 5. Anotemos la variable хз а la 
izquierda de la igualdad enfrente del 5 suprimiendo el menos que 
se encuentra en el cruce de esta línea y la columna 26. Suprimimos 
también el menos del cruce de la segunda línea y la columna 26, 
y de la segunda linea que corresponde al número 7 restamos 
la línea correspondiente al múmero 5. En la segunda línea queda 
el número 2. Eso significa que del segundo punto de partida el 
producto no está todavia distribuido entre los puntos de destino. 
Por eso, en la línea de los valores de (у de la segunda zona 
hallamos el número siguiente por su valor. Ese es el número 7 
que se encuentra en la columna de la variable хуу. Revisando 
la columna encontramos el menor de los números de la izquierda 
que se encuentran al nivel de los menos de la columna (20, 2). 
Suprimimos los menos en la columna 21, anotamos la variable xz, 
a la izquierda de la igualdad: enfrente del 2 y restamos de la 
séptima línea la segunda. Con esto acaba la distribución de las 
unidades de la segunda línea. Continuamos este proceso hasta 
la distribución de las unidades de la sexta línea. Como resultado 
obtendremos la primera solución que coincide con la solución 
anotada en la tabla 24. La tabla grande se ha transformado 
ahora en la tabla 25. 

De la justeza de la solución nos convencemos por la falta 
de los "+" y los “—" en una de las líneas. En esta linea, 
en el caso de que no haya errores en la distribución, tiene que 
haber sólo ceros. En nuestro ejemplo, la última línea ha resultado 


ser de este tipo. En virtud de la condición Ў a, = У. bj, una де 
a 


las ecuaciones del sistema (27) debe ser consecuencia de las demás. 

Las columnas que corresponden a las variables básicas en la 
tabla están señaladas con un ráyado vertical y las columnas en las 
que su ty > 28, con un rayado inclinado. Estas últimas se excluyen 
de la revisión posterior. En la tabla 25 están indicados los modos 
de composición de una mejor variante en comparación con la 
primera. Realmente, nosotros deseamos hallar una solución tal en 
la que la variable que se encuentre en la misma columna que 
ty = 28, es decir la xg, sea igual а cero. Ahora bien, el contenido 
de la sexta línea de esta tabla muestra que la disminución del 
valor de la variable ху puede ser obtenida a cuenta del aumento 
de хеј о de xes. Por cuanto a la variable xg, le corresponde 
tor = 15 y a la variable хез, tes = 23 haremos la disminución 
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de xg a costa del aumento de xg. Revisando la columna 61 
vemos que el menor número que se encuentra en la parte izquierda 
de la tabla enfrente del “=” de esta columna es el valor xg = 5. 

Ahora realicemos la transformación de la última tabla de la 
siguiente manera: excluyamos las columnas que tienen el rayado 
inclinado y la columna que corresponde a la variable xs, (о bien, 
lo que es lo mismo, al número пу, igual a 28). En logar de xo 
escribamos хе. Siguiendo, de las líneas que tienen signo menos 
en la columna 61 restamos la línea que corresponde a la variable 
хвл, y a las líneas que tienen signo más en la columna 61, les 
sumamos la línea que corresponde a esa misma variable. Después 


de esto desaparecerán todos los “—" y los “+” de la columna 61. 
Tabla 26 
RA aa] 
сыга = EE 

a3 = AAA =] 
иш = GOOR 7] 
1 E 
FE E SH +] 

Eje al 
+ 5 == === Û 
zise] O Ш 02 000580 


Como resultado obtendremos la tabla 26. De esta tabla se 
deduce que ya hemos compuesto la segunda solución para la 
realización de la cual se necesita bastante menos tiempo: en lugar 
de 28 horas, sólo 21. (Observemos que si hubiésemos comenzado 
la reducción de la variable xs, a cero no a costa del aumento 
de ха, sino del aumento de xes, hubiese sido compuesta una 
variante de solución con un tiempo máximo = гез = 23. Después, 
tachando todas las columnas соп (у > 23, hubiésemos tenido que 
hacer otro paso para componer una tercera solución en la que 
el máximo de tiempo fuese t = 21). 

Enunciemos las variables, el cambio de cuyo valor ha llevado 
a la segunda solución: хез, xs, X47, Хуу, Xas, Xer. Estas variables 
han formado una cadena cerrada (tabla 27). Numeremos las celdas 
con las variables que forman la cadena tomando como primera 
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la celda con la variable хет (la que se tiene que excluir). Esta 
cadena tiene la propiedad de que en sus celdas impares infalible- 
mente se encuentran elementos x-seleccionados у en las pares, los 
valores гу < tj“, donde г es el valor de ty que debe de ser 
excluido (en este caso (2 = 28). Llamemos aligeradora a esta 
cadena. Después de hallar la cadena aligeradora para formar una, 
nueva solución traspasamos por ella una cantidad de unidades 
de carga que sea igual a la minima entre ellas en la semicadena 
impar, En la tabla pequeña esta composición nos lleva a una 
segunda solución que coincide con la segunda solución compuesta 
en la tabla grande. Excluyamos de la investigación posterior 1 
columnas (celdas) con valores ty > 21. Asi obtendremos la tabla 


Tabla 27 


pequeña 28 y la tabla grande 26 (contando que las columnas 41, 
55 y 65 están tachadas). Surge una pregunta: ¿habrá alguna solución 
que se realice en menos de 21 horas? Para dar contestación 
a esta pregunta examinamos la tabla 26. En la línca de los 
valores ty hallamos el mayor de ellos, el 21. A este valor t4 = 21 
le corresponde la variable ху, en la linea de la cual hay sólo 


signos “+”, La falta de los “—" і lidad de 


la posibilidad de componer una cadena aligeradora en la tabla 
pequeña 28. 

Para aclarar mejor el método de cálculo de la variante óptima 
por el criterio del tiempo, hemos empleado tablas grandes 
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y pequeñas. En realidad, al resolver tales problemas hace falta 
sólo una tabla grande o pequeña hecha a lápiz. Los pasos de una 
variante a otra se efectúan en una tabla empleando lápiz y goma 
de borrar. 


Tabla 28 


сатни 

Las reglas para el cálculo de la variante óptima рог el criterio 
del tiempo cuando se emplea una tabla pequeña se reducen 
a lo siguiente: 

1. Anotamos las condiciones del problema en una tabla pequeña. 

2. Encontramos la primera solución (por ejemplo, con el proce- 
dimiento indicado anteriormente). 

3. Determinamos el ган que corresponde a esta solución. 

4, Tachamos еп la tabla todas las celdas que tengan 1 > timis, 

5. Investigamos la posibilidad de componer cadenas aligeradoras 
formadas de шума“ con los elementos que quedan en la tabla. 

6. Si aparecen cadenas aligeradoras componemos una solución 
nueva. 

7. Si no es posible aligerar por completo la celda que 
corresponde а гё (reducir la respectiva xy а cero) entonces 
la solución que contiene цучх es la óptima. 

8. Si se aligera por completo la celda que contiene шуме, 
entonces hallamos {fri en la nueva solución obtenida. 

Así hemos cumplido el primer paso dirigido a la obtención 
de la solución óptima. Si en este caso no se ha llegado a la 
solución óptima entonces hay que hacer un segundo paso, 
cumpliendo de nuevo las reglas anotadas desde el punto 4. 

Para obtener la solución óptima habrá que realizar un número 
de pasos finito. 

Al emplear estas reglas, la única dificultad consiste en la 
composición de las cadenas aligeradoras. Con respecto а esto, 
el empleo de la tabla grande, cuyas reglas de transformación 
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descritas anteriormente con suficiente detalle, da la 
idad de llegar sin grandes impedimentos a la solución 
óptima, 

Se podrían exponer también otros procedimientos de resolución 
de este problema, 

Cuando т y п son pequeños, estos problemas pueden ser 
resueltos con papel y lápiz Si т y n son grandes entonces es 
necesario, lo mismo que en los problemas que se resuelven por 
el criterio del costo, el empleo de las máquinas computadoras 
electrónicas. En este caso el tiempo de máquina de resolución 
de estos problemas, para los; mismos valores de т y п, coincide 
aproximadamente con el tiempo de resolución de problemas por 
el criterio del costo. 


$ 17. RESOLUCIÓN DE LOS PROBLEMAS 
DE TRANSPORTE TOMANDO EN CUENTA 
EL TIEMPO Y EL COSTO 


En la práctica se presentan condiciones en las que es más 
conveniente la realización de un plan de transporte que sea un 
promedio entre el óptimo por tiempo y el óptimo por costo. 
Supongamos, para hacer más sencilla la descripción, que el costo 
es proporcional al tiempo de transporte, o sea су = Ки. 

En la tabla 29 se presenta la variante óptima, por el criterio 
de costos, del problema visto anteriormente que se cumple en 
28 horas. Los gastos de transporte de esta variante son C = 
= К (7.15+6.5+10.2+11.20+20.13+21.12+5.9+ 12-214 
+ 14.12 + 16-11 + 28.5) = К. 1668 unidades de costo. Los gastos 
en el transporte que corresponden а la variante óptima рог el 
criterio del tiempo (tabla 28) son de С = K (7.15 + 7-2 + 6.54 


Tabla 29 
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+ 11-134 20-13 + 21.12 + 21.7 + 5-2 + 12-28 + 14.12 + 
+ 15۰5 + 16-11) = 1716K unidades de costo. Haciendo un gasto 
complementario, igual a 1716K — 1668K = 48К acabaremos las 
operaciones de transporte de mercancias 7 horas antes en 
comparación con el plan de transporte determinado por la variante 
óptima, calculada por el criterio del costo. 

Del ejemplo examinado se deduce que, admitiendo un gasto 
complementario en el transporte relativamente pequeño, se puede 
reducir considerablemente el tiempo de las operaciones de transporte 


Tabla 30 
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de los productos. En los casos en los que se dan los plazos del 
cumplimiento de las operaciones de transporte se puede, empleando 
los métodos descritos, hallar la variante óptima por el criterio 
de costos que se realice en el margen del plazo dado. Si la pérdida 
de tiempo es inadmisible, entonces calculamos la variante óptima 
sólo por el criterio del tiempo. Si la variante obtenida con tgtno 
es la mejor por los costos, entonces empleando el algoritmo del 
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cálculo de la solución óptima por el mínimo de la forma lineal, 
podremos obtener la óptima por los costos entre todas las óptimas 
por el tiempo (o sea, entre todas las que se realizan en un mismo 
plazo, igual а гун}. Como se че en el ejemplo examinado, la 
solución óptima por el tiempo presentada en la. tabla 30 по es 
óptima por los costos. Efectivamente, reduzcamos a cero los 
elementos x-seleccionados (al reducir а cero los elementos 
x-seleccionados, по зе toma en cuenta el valor de los elementos 


Tabla 32 
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que se encuentran en las celdas tachadas, tabla 31) y como resultado 
obtendremos una tabla en la que hay un número negativo (— 14). 
Transformando otra vez la tabla 31 obtendremos la tabla 32. 
En la tabla 32 ya transformada los elementos x-seleccionados son 
iguales a O y los demás son no negativos. Por eso, la variante 
de solución compuesta es la óptima por el criterio del costo. 
El costo de la realización de esta variante, que se lleva a cabo 
también en 21 horas, es de С = 1688K lo que es mucho menos 
que 1716К. 


En los últimos años, los métodos de programación lineal 
encuentran cada vez mayor aplicación en la economía, la técnica, 
el arte militar, y en otras esferas. La teoría y los métodos de 
programación lineal se perfeccionan constantemente y permiten 
solucionar otros nuevos problemas. El desarrollo acelerado de la 
técnica de computación ha hecho posible la solución práctica de 
cualquier problema de programación lineal. 

La continuación del desarrollo de los métodos de programación 
lineal y su empleo para la resolución de problemas de la economía 
nacional va a favorecer la mejora de la organización y de la 
planificación de la industria nacional. 
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